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ABSTRACT 


Automatic flare and decrab control laws were developed for NASA's experi* 
mental Twin Otter. This light wing loading STOL aircraft was equipped with 
d1 'ect lift control wing spoilers to enhance flight path control. Automatic 
I tnding control laws that made use of the spoilers were developed, evaluated 
In a simulation and the results conq)ared with these obtained for configurations 
that did not use OLC. The spoilers produced a significant Improvement In per- 
formance. A simulation that could be operated faster than real time In order to 
provide statistical landing data for a large number of landings over a wide 
spectrum of disturbances In a short time was constructed and used In the evaluation 
and refinement of control law configurations. A longitudinal control law that 
had been previously developed and evaluated In flight was also simulated and Its 
performance compared with that of the control laws developed In this study. Run- 
way alignment control laws were also defined, evaluated and refined to result In 
a final recommended configuration. Good landing performance, compatible with 
Category III operation into STOL runways, was obtained. 
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1.0 INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 


The Ames research center of NASA conducted a series of investigations 
to generate and verify through ground based simulation and flight research 
a di.ta base to aid In the design and certification of advanced short takeoff 
and landing (STOL) aircraft. A modified de Havllland Canada DHC-0, Twin Otter, 
with direct lift wing spoilers was evaluated as a representative light wing 
loading STOL aircraft. The Augmentor Wing Jet STOL research airplane and the 
Quiet Short-Haul Research Aircraft (QSRA) are representatives of the propulsive 
lift technology. One portion of this program Is concerned with obtaining techni- 
cal Information on automatic landing systems for STOL aircraft Including flight 
path control performance and touchdown state dispersion In the presence of en- 
vironmental disturbances. As part of this program, Lear Slegler's Astronics 
Division developed autorriatic landing control laws for these experimental 
aircraft. The results of Lear Slegler's Twin Otter work are presented in this 
report and previous studies are documented as References 1, 2 and 3. 

The technology for the development and certification of Category III auto- 
matic landing systems for conventional takeoff and landing (CTOL) jet transports 
is well developed and documented, as noted In References 4 to 6 for one commer- 
cial aircraft and Reference 7 for the FAA requirements. No comparable technology 
exists for automatic landing systems for STOL airplanes. 

The objective of the automatic landing work reported here is to gain under- 
standing of the problems impacting the design of light wing loading short-haul 
airplanes that are to be landed automatically on STOL runways in adverse weather 
conditions. This understanding was attained by a limited coverage of important 
elements that are normally included In the certification process of a CAT III 
automatic landing system for CTOL airplanes with major emphasis on fault-free 
performance. The control law development concentrated on the final approach to 
touchdown phase of the landing. 

Several longitudinal automatic landing control laws, with and without the 
use of the OLC spoilers, were developed as part of this study and their perfor- 
mance with atmospheric disturbances and landing aid noise was evaluated deter- 
ministically and statistically in a simulation. Another longitudinal control 
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law that had been defined Independently of this study and evaluated In flight 
was evaluated In the simulation for comparison against the newly developed 
control laws. Lateral and directional control laws were also defined during 
this study and evaluated In the simulation. 

1.2 SUWARY 

This report describes the development of a family of automatic landing 
system control laws and shows that this type of control law Is capable of 
meeting requirements like those applied by the FAA to CTOL automatic landing 
systems. The results presented In this report are derived from simulation. 

The report contains seven sections describing the development and evalViatlon 
of the automatic landing control laws. Section 2 Is a brief description of the 
Twin Otter Airplane, the STOL approach conditions, the airplane's controls and 
Its avionic system. Section 3 describes the design and evaluation process em> 
ployed In this program. Section 4 contains a detailed description of the final 
longitudinal and lateral control laws that emerged from this program. Section 5 
describes the longitudinal landing performance results that were obtained In 
the simulation. It provides data for the deterministic and stochastic perfor- 
mance of the control laws that were evaluated. Section 6 describes the simula- 
tion landing perfonnance result* for the lateral /directional axis. The impact 
of various error sources Is analyzed and tradeoff studies leading to the final 
recommended control law are described. The conclusions derived from this work 
are presented in Section 7. Appendix A Is a summary of airframe, controllers, 
sensors and disturbance mathematical models that were used In the simulation. 
Appendix 6 contains backup data for the longitudinal simulation results that 
are presented in Section 5 and Appendix C contains backup data for the lateral/ 
directional results of Section 6. 

2.0 THE RESEARCH AIRPLANE AND THE APPROACH CONDITION 

A three view drawing of the experimental OHC-6, Twin Otter, 1$ shown as 
Figure 2-1. The airplane had been originally designed as a light wing loading 
STOL airplane and therefore xhe only major airframe modification in the experi- 
mental version was the incorporation of the wing spoilers that are shown in 
Figure 2-1. 
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The nominal landing approach condition of the Twin Otter Is compared In 
Table 2-1 with that of the Auqmentor Ming Research Aircraft which Is a oowered 
lift STOL airplane and the Lockheed L-1011 which Is a typical examole of a CTOL 
transport airplane. 

table 2-1 STOL VERSUS CTOL LANDING APPROACH CONDITIONS 



AUG WING 

TWIN OTTER 

L-1011 

A1'^rj.?ed, kt 

70 

70 

135 

Glloe Slope Angle, deg 

7.5 

fi.O 

2.75 

Thru-t Inclination, deg 

90 

SO 

0 

Wing Loading N/m2 fpsf) 

2378 (49.71 

1254 (26.21 

4359 (91.11 

Approach Lift Coefficient 

3.0 

1.5 

1.5 

Approach /Full Thrust, percent 

85 

10 

25 

Lift, Aerodynamic, percent 

40 

100 

100 

Cold Flow, percent 

40 

0 

0 

Hot Thrust, percent 

20 

0 

0 


NOTE: These numbers are representative and aoproxlmate; they are given In order to 
hlqhllaht the differences between the three aircraft rather than to provide 
exact data for each one. 

The slow and steeo approach Is common to both STOL aircraft and It Is In 
contrast with the fast and shallow aooroach of the CTOL airplane. The llaht wina 
loading Twin Otter, like the L-1011, derives essentially all its lift from 
aerodynamic sources, flies at a conventional lift coefficient value and with a low 
power setting. The relatively high thrust Inclination of the Twin Otter is a 
result of the propeller flow acting on the wing and flap. Unlike the Au<inentor 
Wing, the Twin Otter operates on the front side of the power curve at 1.3 times the 
stall speed and therefore the conventional control technique of using the elevator 
for flight path control and throttle for speed Is effective for this airplane. 
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The spoilers shoi<n In Figure 2-1 are the only non standard control on the 
Twin Otter. Each wing has a total of six panels, three on the upper surface and 
three on the lower one. The four Inboard panels modulate lift and drag and the two 
outboard spoilers enhance roll control hut were not used In this study. Of the 
four lift/ drag spoilers, only the outboard uoper panel was used since It was the 
most effective one. Flight tests show that the Inner upper panel Interacts with 
the horizontal tall and causes excesslye buffetting and the lower panels haye 
little aerodynamic effect. When used, the soollers are biased at 20® and the 
airspeed ^s Increased by four knots fto 75 knots at nominal weight) to compensate 
for the attendant Increase In stall speed. The spoilers can he modulated between 
zero and 40* In the glide slope track phase, proyiding a direct lift control (DLC) 
authority of ±0.13 g. A buffetting tendency Is associated with spoiler deflections 
of more than 40* but during the flare, deflections up to ^n® are allowed. 

Lateral/d1?*ect1onal control Is conyentlonal with ailerons for roll and rudder 
for yaw. 

It'e airplane Is eoulpoed with the STOLANO digital aylonics system (Reference 
8) providing versatile navigation, guidance, control and display functions. 

During the first test period, two successive microwave landing systems were 
used for approach guidance, providing azimuth, elevation and distance Information. 
The first of these was an experimental system called MOOILS which had narrow 
azimuth coverage and a relatively granular azimuth signal. The second was a 
prototype of the microwave landing system that has been adapted by the Federal 
Avi ati on Adml ni strati on . 

3.0 DESIGN AND EVALUATION fCTHOOS 

3.1 SIMULATION 


The development of a simulation allowing the collection of automatic landina 
performance with deterministic and stochastic disturbances Is an Important part of 
the certification process of a Category III automatic landing system for .. CTOL 
aircraft, as described In References 1 through 6. Such a simulation of the Twin 
Otter was developed for thl- study as a major design and evaluation tool for 
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the control laws that were developed as part of this study. This simulation was 
also used for the evaluation of longitudinal control laws that were developed 
independently of this study. The best available airframe data were used in the 
simulation but a rigorous correlation between the simulation and flight test 
results was not included in this study. Also, the control laws that were developed 
during this study were not evaluated in flight as a result of flight time 
constraints. Thus, all landing performance results that are given in this report 
were obtained from the simulation. 

The simulation that was used as the major tool for synthesizing and 

evaluating the automatic landing control laws was designed to operate ten times 

faster than real time in order to allow the collection of statistical data. 

Mathematical models of the airframe, controllers, sensors and the environment were 

assembled and used in the simulation. The normal set of uncoupled, linearized, 

small perturbation equations of motion were used in separate longitudinal and 

lateral simulations. Longitudinal dynimics were included in the lateral simulation 

to the extent necessary to account for the ground speeds associated with different 

headwinds. Lift, pitching moment, and drag variation due to ground effects were 

also included. Controller dynamics were modeled, including rate and position 

limits. Sensor dynamics and error models which contribute to landing dispersions 

were also induced, such as radar altimeter dynamics and offsets, and dynamic and 

static vertical gyro and accelerometer errors. MLS noise was modeled and included 

in the simulation. Winds, shears and lurbulence consistent with the definitions in 

* 

the FAA Advisory Circular 20-57A (Reference 7) were used. 

For statistical data collection, the simulation was run in fast time 
repetitive operation mode, starting at 152.4 m (500 ft) above the runway with the 
airplane stabilized on the glide slope or localizer, and terminating at touchdown. 
The 30.5 meters (100 foot) approach window states were recorded, as were the 
touchdown states: vertical and lateral velocity, touchdown point on the runway, 

and pitch, roll and heading angles. Data were taken for various levels of 
environmental disturbances, ar'! system errors, covering a wider range than possible 
in flight. Probability distributions were generated for all touchdown state 
variables. 

Appendix A contains a detailed description of the models used in the 
simulation. 
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3.2 THE ^JOl RUNWAY 


The Twin Otter was flown by NASA Ames Research Center at Crows Landing Naval 
Auxiliary Landing Field (NALF) in California. The flight test landings were msde 
on a simulated 518.3 by 30.5 meters (1700 by 100 feet) ST(M. runway with boundaries 
oainted, in accordance with Reference 9, on a longer and wider runway. The runway 
geometry is shown in Figure 3-1. The 5® Glide Path Intercept Point is 71.6 meters 
(235 ft) beyond the threshold. All landing distance results in this report are 
referenced to the SPIP. 


The marked touchdown ?one extends from 19.8 m (65 ft) to 80.8 m (265 ft) 
bevond the SPIP. 



Figure 3-1. STOL Runway Geometry 
3.3 LANDING PERFORMANCE OVERVIEW 

A one signa landing dispersion summary, comparing results obtained from 
simulation with limiting shear winds, moderate turbulence and beam noise for the 
Twin Otter with comparable Augmentor Wing results, is given in Table 3-1. The two 
aircraft differ significantly in their method of control in the loioitudinal plane 
but are similarly controlled in the lateral /directional axis. The landing 
performance results are guite similar, however, in both axes. 

Limiting winds refer to the limit values considered in the atmospheric 
disturbance model described in Appendix A; i.e., 25 kt headwinds and 10 kt 
tailwinds together with corresponding values of turbulence and shear. These levels 
should result in rather conservative estimates of aircraft dispersion due to this 
source but they provide a useful basis of comparison of performance between 
different aircraft and/or control systems. 



TABLE 3-1 ONE SIGMA LANDING DISPERSION SUMMARY 


TWIN OHER 


NO OLC 


h m/sec (fps) 0.213 (0.70) 0.320 (1.05) 


X m (ft) 


30.5 (100) 36.6 (120) 


0 deg 


y m (ft) 


1.49 (4.9) 


y m/sec (fps) 


0.26 (0.85) 


deg 


AUG WING 


NO OLC 


0.229 (0.75) 0.381 (1.25) 


32.0 (105) 47.2 (155) 


1.28 (4.2) 


0.15 (0.50) 




NOTES: 

1) The results shown are the dispersion between the mean and the one sigma point. 

2) Simulation data with 70% limiting shearing headwind, 30% limiting shearing tail- 
wind, the corresponding moderate turbulence, and beam noise. 

3) Twin Otter longitudinal results are for the constant flare height low gain 
configuration, as defined in Section 5. Three and Two Control Augmentor wing 
results are given, as described in Reference 1. 


3-4 



Sink rate and range control for both aircraft benefit from the use of 
direct lift control. The Augmentor Wing no DLC longitudinal performance is 
Z0% to 30X worse than the comparable Twin Otter results. Without OLC, flight 
path angle is controlled with thrust for the high wing loading Augmentor Wing 
airplane, and with attitude for the light wing loading Twin Otter. The results 
suggest that the latter Is more effective for sink rate and range control. 

With OLC, both aircraft have similar performance as their control bandwidths 
are equalized. Pitch attitude dispersion is, somewhat surprisingly, similar 
for the two airplanes. Lateral results are also similar, with the Augmentor 
Wing being somewhat better In lateral displacement and velocity whereas the 
Twin Otter is aligned more accurately with the runway heading with somewhat 
less roll activity. 

The Twin Otter landing performance results are described In detail in 
Section 5 for the longitudinal axis and In Section 6 for the lateral. Aug- 
mentor Wing results are given in Reference 1. 
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4.0 DESCRIPTION OF CONTROL LAWS 


Two sets of control laws for automatic glide slooe track and flare were 
evniuated in this study. One set was defined as oart of the 'ork reoorted here and 
.c utilizes a constant flare initiation height and closed looo oath control about a 
oredetermi ned oitch trajectory command. The other set had been defined orior to 
this study and it uses a variable flare height and oredictive oitch and elevator 
commands in addition to the closed loop control. The variable flare height control 
laws »ere evaluated in flight and in simulation whereas the constant flare height 
mechanization was evaluated in simulation only. Both sets of control laws were 
evaluated with and without the use of the OLC sooilers. Both control laws are 
described in Section 4.1. Evaluation results are given in Section 5. 

Control laws for automatic localizer track and runway alignment were also 
defined and evaluated in simulation as oart of this study. The final configuration 
that emerged from tradeoff studies is described here in Section 4.2. Tradeoff 
studies and evaluation results are oresented in Section f. 

4.1 GLIDE SLOPE TRACK AND FLARE 

4.1.1 CONSTANT FLARE HEIGHT • 


The constant flare height longitudinal control laws that have been developed 
for the Twin Otter as oart of thi^ study are shown in the block diagram of Fioure 
4-1 and the numerical values of gains an constants are defined in Table 4-i. A 
conventional, front side of the power curve, control technioue is used on this 
light wina loading STOL airplane. The elevator is used for the control of flight 
oath angle as well as attitude stabilization and control. The throttles are used 
for speed control. The DLC spoilers, when used, assist in the control of flight 
oath angle. 

Four filters are used to smooth inputs into the control law in order to 
reduce noise induced control activity. The block diagrams of the filters are given 
in Figure 4-2. Raw airspeed is blended with longitudinal acceleration in a 
complementary filter to produce a smoothed estimate of airspeed which is used to 
drive the throttles. Raw glide slooe error, computed from elevation and range 
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RADAR ALTIMETER 











TABLE 4-1 LONGITUDINAL GAINS AND CONSTANTS 
CONSTANT FLARE HEIGHT 


GAINS WITH SPOILERS WITHOUT SPOILERS 



deg/deg/ sec 

1.0 

1.0 


s 

deg/deg 

2.0 

2.0 


hi 

deg/sec/deg 

0.50 

0.50 

S 

deg/m/sec (deg/fps) 

3.28 (1.0) 

3.28 

(1.0) 

S‘e 

deg/m/sec^ (deg/fps^) 

0.656 (0.2) 

0.656 

(0.2) 

Sle 

dec/ sec # deg/ sec « 
m/sec ^ rps ' 

0.656 (0.2) 

0.656 

(0.2) 

S 

m/sec/m (fps/ft) 

0.50 

0.50 

Si 

m/sec^/m (fps^/ft) 

0.050 

0.050 

s 

cm/kt (1n/kt) 

0.323 (0.127) 

0.323 

(0.127) 

Si 

cm/sec/kt (In/sec/kt) 

0.0323 (0.0127) 

0.0323 

(0.0127) 

St 

cm/m/sec (in/fps) 

2.50 (0.30) 

2.50 

(0.30) 

Sp 

deg/m/sec (deg/fps) 

22.97 (7.0) 

0 


Ss 

(fps/fps2, 

0.50 

- 


So 

- 

0.20 

0 


constants 





^GS 

Q 

(ft) 

30.48 

(100) 

30.48 

(100) 

^FL 

m 

(ft) 

12.19 

(40) 

18.29 

(60) 

*'SFL 

m 

(ft) 

15.24 

(50) 

15.24 

(50) 

^FL 

m 

(ft) 

12.19 

(40) 

12.19 

(40) 

N’CF 

m 

(ft) 

12.19 

(40) 

9.14 

(30) 

^0 

m/sec (fps) 

0.853 

(2.8) 

1.219 

(4.0) 

. c 
^0 

deg 


4.5 

3. 

5 


kt 

- 

14.8 

-14. 

8 


TIME CONST^XNTS 


^awo 

■^SWO 

TSP 


5&C 

%sc 

sec 


sec 


A0 


LIM 


15.0 

0.10 

5.0 


0.5 


*10* 9 h > h 
* 10 * 


- "epL 
h 




15.0 

0.10 
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'information, is combined with washed out normal acceleration in anoth<;r comple- 

f 

mentary filter to produce an estimate of deviation from the glide-slope. A 
simple 0.1 second first order lag is used to filter the output of the radio 
altimeter and to derive sink rate. The derived sink rate is blended with 
washed out normal acceleration and the resultant sink rate signal is used in 
the flare. 

Pitch attitude and rate feedbacks to the alevator are used for attitude 
stabilization and control (Tigure 4-1). A pitch error integrator is used to 
maintain elevator trim. Soioothed glide-slope error is summed with the output 
of the glide-slope Integrator to produce a sink rate error v'hich is sunned 
with washed out and lagged normal acceleration and commands attitude corrections. 
Washed out normal acceleration and raw glide slope error are used as inputs to 
the glide slope integrator such that it provides glide slope damping and helps 
to null the steady state error. While in the glide slope track mode, the radio 
altimeter based sink rate signal is subtracted ahead of the flare scheduler and 
added past it. Since this scheduler is at unity gain prior to flare the net 
result is that this path has no effect on the system. Attitude reference is 
generated by the sink rate error integrator while in the glide slope track 
mode. A transition from the glide slope track to the flare mode occurs be- 
tween a gear height of (30.48 m or 100 ft) and h^p^ (15.24 m or 50 ft). 

The glide slope error signal is faded out and the airplane maintains the glide 
slope sink rate. At hgp|^ an attitude change is conmanded linearly with de- 
creasing altitude from the approach to the touchdown value, as shown in Figure 
4-3. The rotation arrests the sink rate and puts the airplane in a proper touch- 
down attitude. Through the flare, derived sink-rate is transitioned linearly 
with decreasing altitude from glide slope to radar altimeter based Information, 
minimizing the impact of terrain irregularities. A straight line ii/fi profile from 
the existing pre- flare sink-rate to the desired touchdown value is commanded in 
the flare as shown in Figure 4-3. This results in an exponential flare, the 
time constant of which is proportional to the slope of the h/li line. The flare 
height is constant at 12.19 m (40 ft). The pre-flare sink-rate varies with the 
wind conditions, resulting in a flare time that is shorter with tailwind than 
with headwind. This variation in flare time tends to compensate for wind in- 
duced touchdown position dispersion. Without spoilers, a constant touchdown 
sink rate coninand of 1.219 m/sec (4.0 fps) is used. With spoilers, however, 
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the flare scheduler qain at touchdown Is non zero 0.2, as shown 1n Table 

4-1) and consequently an Incremented touchdown sink rate command, that Is a 
function of pre-flare sink rate. Is Introduced. The sink rate error Is used to 
modulate the linear pitch attitude command that Is shown In Flqure 4-3. The 
allowed nose down command Is reduced from -10® at pitch flare helqht to zero at 
touchdown, as Indicated by 1n Figure 4-1 and Table 4-1. 

The throttles are used to maintain the approach airspeed while trackinq the 
glldeslooe. A speed reduction command Is applied during the flare. Increasing 
linearly with decreasing height. Sink rate error Is cross fed to the throttles 
during the flare In order to enhance vertical control and compensate for the 
limited ability to use nose down commands close to the ground. 

When the DLC spoilers are used they are driven with a blend of sink rate 
error and washed out and lagged normal acceleration. The spoilers are used for 
short term flight oath angle corrections both In glide slope track and flare. The 
spoilers allow a faster modulation of normal acceleration th^ pitch changes do and 
therefore they allow an increase of system bandwidth and Improve the rejection of 
disturbances. In glide slope track, the spoiler command Is belno washed out with a 
five second time constant In order to preserve their full dynamic range to counter 
rapid disturbances while handling slow trim changes with pitch attitude. In the 
flare, the spoiler washout Is eliminated such that the full available spoiler lift 
modulation capability can be used and held If needed. 

4.1.2 VARIABLE FLARE HEIGHT 


The variable flare helqht control law had been designed prior to this study, 
excluding the spoiler control. It has been tested and refined In flight. One task 
of this study called for simulating this control law in order to evaluate its 
performance in comparison with the constant flare height control law. Also, the 
spoiler control law that had been developed In this study for the constant flare 
height mechanization was adapted and tested In flight with the variable flare 
helqht svstem. Elevator and throttle control of the variable flare helaht law are 
defined bv the block diagram '•f Flqure 4-4. The spoiler control for this svstem is 
shown In Figure 4-^ and the numerical values '^f the associated oalns and constants 
are defined In Table 4-1 1. 
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Variable Flare Height 

























Control Block Diagram ~ Variable Flare Height 











The allocation of controllers is the same as in the previously described 
constant flare height control law The elevator is used for attitude stabilization 
and control and for flight path angle control. The throttles are used to maintain 
airspeed while tracking the glide slope. The spoilers assist in controlling flight 
oatii angle when used. Pitch attitude and rate are tsed as feedbacks to the 
elevator for oitch stabilization. The elevator gain is scheduled with dynamic 
pressure and the elevator and pitch rate gains are higher in the flare than in 
glide slooe track, as can be seen in Figure 4-4. Smoothed glide slope error, its 
integral and estimated glide slope rate generate pitch attitude commands to track 
the glide s1- .n’. At flare, the last glide slooe track oitch command is held and 
retained. F ;.e attitude commands come from two sources. , closed loop command, 
®c, controlling sink rate as a function of gear height and a predictive 
command, 9p^ that is oroportional to the difference between the commanded 
touchdovm attitude and the approach attitude and is increased linearly with 
decreasing altitude. A direct predictive elevator command, is also 
computed based on the difference between the commanded touchdown attitude and the 
pre-flare attitude, pre-flare airspeed and ground speed, as shown in Figure 4-4. 

PI are height is computed as a function of sink rate as shown in Figure 4-4. 

An increase in flare height with tailwind is computed using the ratio of airspeed 
to groundspeed. Flare height as a function of sink rate is given in Figure 4-6. 

One line is shown for headwind or zero wind and another one for a 10 knot steady 
tailwind. The closed loop commanded flare trajectories for a 25 knot steady 
headwind, zero headwind, and a 10 knot steady tailwind are also shown in Figure 
4-6. These trajectories are the loci of o' » o. Above the trajectory a nose 
down command is computed and below the trajectory nose up is commanded. The 
commanded touchdown sink rate is constant at 0.686 m/sec (2.25 fos). The nominal 
'■“e-flare sink rate for an approach speed of 71 knots and a glide slooe of 6* is 
also shown in Figure 4-6 for a 25 knot steady headwind, zero headwind and 10 knots 
steady tailwind. In all cases flare height is above the commanded closed loop 
trajectory. This is done in order to allow time to "turn the corner" and minimize 
the undershoot of the sink rate trajectory. The initial rotation is provided bv 
the pitch and elevator predict terms while e' is inhibited during the initial 
1.5 seconds of the flare if nose down is being commanded, as shown in Figure 4-4. 

A transition from the flare height line to the closed loop control line is computed 
through use of the "flare switch" circuit shown in Figure 4-4. 9' is faded in 

at flare initiation. The gain from 0' to oitch attitude is oroportional to 
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the ratio between pre-flare ground speed to airspeed. Below 3.81 m (12.5 ft) 
higher gain is used for nose up conmands than for nose down. Vertical acceleration 
feedback is added to the pitch conmand during flare to improve damping. 

During glide slope track the throttles are used to control airsoeed. An open 
loop throttle trim position Is supplied from trim tables. Filtered airsoeed error, 
the integrol of raw airspeed error and longitudinal acceleration are summed to 
provide a closed loop throttle connand. At flare, this conmand is held constant 
and a constant rate retard is computed from pre-flare sink rate and airspeed. The 
fact that the throttles are not used actively through the flare is considered as a 
disadvantage of this mechanization in comparison with the constant flare height 
control law. 

When the spoilers are used they are biased to 20*. As in the constant flare 
height control law, they are used to provide rapid short term flight oath angle 
corrections. While tracking the glide slope, filtered glide slope error, estimated 
glide slope deviation rate and washed out vertical acceleration are summed with the 
appropriate gains to generate the spoiler conmand (Figure 4-5). This conmand is 
washed out with a 3 second time constant in order to maintain the full steady state 
lift modulation capability. At flare, the glide slope conmand is faded out and the 

flare conmand, made of 9*^ washed out vertical acceleration, is faded in. 

The flare command to the spoilers is not washed out. 

4.2 LOCALIZER TRACK AND RUNWAY ALIGNMENT 


The sensitivity to external disturbances and large crab angles typical of 
STOL aircraft places stringent demands on the design of the automatic landing 
system. The accurate control of aircraft position and heading required Just prior 
to touchdown makes the runway alignment the most critical phase of the lateral 
landing control problem. Several candidate control laws were developed as oart of 
this study and evaluated by simulation. All the control laws developed in this 
study used toward slip for runway alignment because this technique was determined 
to be superior in a study that had been done for the Augmentor Wing Aircraft 
(Reference 3). Tradeoffs were made in other areas and the results are presented in 
Section 6. The control law that produced the best results is described in this 
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section. Figure 4-7 Is a block diagram of the localizer track and runwav alignment 
control law. The numerical values of the associated gains and constants are 
defined In Table 4-1 II. Roll control on the Twin Otter alrolane Is conventional 
with the control wheel mechanically linked to the ailerons. The lateral control 
law outout commands a wheel position for roll control. Raw localizer lateral 
displacement computed from azimuth angle deviation and range, is blended with cross 
trick acceleration In a complementary filter. The estimated localizer deviation 
and Its rate are used to command bank angle. Additional pseudo rate is provided bv 
filtering cross track acceleration with a 20 second first order lag. This helps 
reduce MLS beam noise Induced control activity. The yaw rate, lateral acceleration 
and bank angle command signals are fed through gains and summed to drive the rudder 
for yaw stability augmentation and turn coordination. Yaw rate and rudder position 
command cross fed to the wheel are used to minimize the degradation of localizer 
tracking accuracy due to the airplane* - nounced dihedral effect. 

A forward slip maneuver Is used for runwav alignment. Reginninq at an 
altitude of 45.7 m fl50 ft), an align command Is switched Into the vaw axis. This 
reference heading command Is reduced from the heading error existing at alignment 
Initiation to zero at 15.2 m (50 ft), yielding an alignment rate which Is a 
function of both initial heading error and aircraft sink-rate. The error from the 
commanded heading trajectory Is Integrated to mai tain the steady rudder reoulred 
during alignment. In the roll axis, the beam computations are maintained to guide 
the vehicle along the desired horirontal oath, with increased oseudo cross track 
rate gain for better control. Wing down compensation Is provided by Inserting 
C'^oss track acceleration and bank angle through a one second lag and a ±5* limit. 
This results in a one second washout on the roll attitude feedback within the ±5® 
bank limit which allows roll attitude to reach the value needed to null the cross 
track acceleration. The 5® alignment limit Is ample for steady crosswind levels 
higher than the 15 knot specifed by the FAA as the limiting case (Reference 7), and 
minimizes the potential for laroe touchdown bank angles. If the limit is exceeded, 
the vehicle will maintain a small crab angle. 
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Localizer Track and Runway Alignment Block Diagram 















TABLE 4- I II LATERAL GAINS AND CONSTANTS 


GAINS 


H 

deg/deg/sec 


0.25 


deg/deg 


4.0 


deg/sec/deg 


0.1 

•^YA 

deg/m/sec^ (deg/fps^) 

5.84 

(1.78) 

S“l 

deg/m/sec^ (deg/fps^) 

1.23 

(0.375) 

SIa 

deg/m/ sec^ (deg/fps^) 

1.23 

(0.375) 

Sl 

deg/m/sec (deg/fps) 

0.820 

(0.25) 

Sl 

deg/m { deg/ft) 

0.197 

(0.06) 

•S.IM 

deg/sec/deg 


10.0 

Si 

deg/sec/m (deg/sec/ft) 

0.0033 (0.001) 

•^rCF 

deg/deg/sec 


4.0 

ScF 

deg/deg 


0.5 

Sy 

deg/m/sec^ (deg/fps^) 

5.38 

(1.64) 


deg/deg/sec 


2.14 

''n 

deg/deg/sec 


4.0 

Sl 

deg/deg 


0.75 

SiL 

deg/ sec/deg 


0.10 

CONSTANTS 



^LN 

m (ft) 

45.72 

(150) 


m (ft) 

15.24 

(50) 

tALIM 

deg 

' 

±5.0 

\lM 

m/sec^ fps^ 

±0.975 

(±3.2) 

*CLIM 

deg/sec 


±3.0 

*CLIM 

deg 


±5.0 

^LIM 

deg 


±2.0 

\lM 

deg/sec 


±4.0 

*ANT 

m (ft) 

(5.97) 

(19.6) 

*■300 

m (ft) 

0.706 

(2.31) 

TIME CONSTANTS 



"1 

sec 


50.0 

^2 

sec 


1.0 

'3 

sec 


20.0 

"4 

sec 


0.10 


sec 


5.0 
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5.0 LONGITUDINAL PERFORMANCE 


The major thrust of this study was to develop automatic landing control laws 
for the Twin Otter STOL airplane and to evaluate the performance of these control 
laws In a fast time simulation that allows the collection of statistical data based 
on a large number of landings. The performance results of the simulation were used 
to refine the control laws In order to optimize their performance. This section 
describes the results obtained from this fast time simulation. 

The constant flare height control law, described In Section 4.1.1 was 
developed and evaluated as part of this study with and without the use of DLC 
spoilers. The variable flare height control law, described In Section 4.1.2 had 
been developed Independently of the work reported here and it has been tested In 
flight and refined by NASA. This control law has been implemented in the fast time 
simulation and its performance was evaluated and compared with that of the constant 
flare height Taw. The results of the simulation evaluation are described in the 
following sections. 

5.1 EVALUATED CONFIGURATIONS 


A total of six glide slope track and flare control !«/ configurations have 
been implemented In the simulation and evaluated during this study as shown in 
Table 5-1. 


TABLE 5-1 EVALUATED CONFIGURATIONS 


FLARE HEIGHT GAINS SPOILERS FLOWN 


1) 

Variable 

Low 

No 

Yes 

2) 

Variable 

Low 

Yes 

Yes 

3) 

Constant 

High 

No 

No 

4) 

Constant 

High 

Yes 

No 

5) 

Constant 

Low 

No 

No 

6) 

Constant 

Low 

Yes 

No 
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Initially, only two constant flare height configurations were developed (Items 
3, 4 In the Table). The variable flare height configurations, however, used 
lower gains in some paths and therefore the low gain constant flare height 
configurations (Items 5 and 6 In the Table) were defined and evaluated In order 
to allow to differentiate between performance benefits that result from control 
law structure and these that are derived from higher gains. The high gain con- 
stant flare height system may be somewhat unrealistic since it has not been 
flown and the simulation math models are probably more benign than the real 
airplane. The low gain constant flare height configuration, however, uses 
gains that are equal or lower than these used In the variable flare height con- 
figuration (which was tested in flight) and therefore it Is very likely to be 
realizable. The high gain constant flare height configuration Is defined In 
Section 4.1.1. It's block diagrams are given In Figures 4-1 and ’4-2. Table 5-II 
defines the numerical values of the gains and constants of the high gain and the 
low gain constant flare height control laws. The time constants, as given in 
Table 4-1, are conmon to both variants. The high gain constant flare height 
system uses the same gains for glide slope track and for flare. The variable 
flare height system uses higher gains for flare In some paths. Therefore, In 
the low gain constant flare height configuration gains are switched in order 
to match the variable flare height values. Table 5-III compares gains of all 
six configurations in all the paths. The gain of pitch rate to elevator in the 
high gain variant of the constaot flare height control law is higher than the 
gain used by the variable flare height law in glide slope track but lower In 
the flare. In the low gain variant, pitch rate to elevator gains are adjusted 
to match these of the variable flare height system. Pitch attitude to elevator 
gain in the nigh gain variant is higher than used by the variable flare height 
law in glide slope track or flare. Again, the gains were adjusted in the low 
gain variant to match these of the variable flare height configuration, as can 
be seen from Table 5-I1I. The high gain constant flare height system includes 
a path of integrated pitch error to the elevator. The variable flare height 
configuration does not have such a path. Thus, this path has been eliminated 
from the glide slope track mode of the low gain variant. However, it has been 
retained in flare as it is needed to compensate for ground effect induced trim 
changes, a function that is handled by the elevator predict terms in the varia- 
ble flare height configuration. Vertical acceleration feedback to attitude 
comnand is used by the variable flare height control law only in the flare. 
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TABLE 5-1 lA LIST OF GAINS AND CONSTANTS - CONSTANT FLARE HEIGHT 

(METRIC UNITS) 





IS. 24 15.24 





- 14.8 - 14.8 



TABLE 5-!IIA GAIN COMPARISON (METRIC UNITS) 
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The high gain variant of the constant flare height configuration uses it both in 
glide slope track and flare. This feedback has been eliminated from the glide 
slope track mode of the low gain variant in order to match the variable flare 
height implementation. The sink rate error to pitch cofftmand gain in the high gain 
variant is slightly lower than any of the egui valent gains used by the variable 
flare height system without spoilers and it was therefore retained in the low gain 
variant. With spoilers, the variable flare height control law uses significantly 
reduced sink rate error to oitch gain and the low gain constant flare height 
mechanization was modified accordingly, as shown in Table 5-III. The constant 
flare height system has a oath of integrated sink rate error to pitch whereas the 
variable flare height control law does not. This oath was retained in the low gain 
version of the constant flare height system* because this integrator is used to 
store the ore-flare pitch trim value. Gains of glide slooe error and it's integral 
to pitch were reduced in the low gain constant flare height system to match 
variable flare height gains. With spoilers, however, the gains used by the low 
gain configuration in these paths are somewhat lower because the same paths are 
used to drive the sooilers and matching the variable flare height spoiler gains 
resulted in lower glide slope and its integral to pitch gains. Longitudinal 
acceleration feedback to the throttles Is used in the variable flare height control 
law but not in the constant fU're height configuration because simulation results 
indicated a degradation in performance If this feedback was included in this 
configuration. Airspeed to throttles gain in Lne constant flare height svstem was 
slightly lower than in the variable flare height configuration to begin with and 
there was therefore no need to reduce It. The gain of the integral of airspeed to 
the throttles was reduced in the low gain configuration. A cross feed of sink rate 
error to the throttles is used hy both variants of the constant flare height 
control law but not by the variable flare height svstem. The svstem that was 
actually tested in flight included software gains in the throttle servo loop that 
oroduced a one second lag from throttle command to actual position. This lag was 
intended to be compatinle with engine 1ooo dynamics. The high gain constant flare 
height control law assumed ..high bandwidth throttle servo (cascaded with the 
engine response as defined in Appendix A). The one second lag in the throttle 
servo was incorporated with the low gain configuration. The gain of vertical 
acceleration to the sooilers is reduced for the flare in the variable flare height 
svstem but the constant flare height configuration retains the glide slooe track 
value in the flare. Sink rate and olid® slooe error gains to the sooilers in the 
constant flare height control law we*‘e close to these used by the variaole flare 
height system and there was therefore no need to reduce them. 
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Table 5-IV summarizes the major differences between the variable flare height 
and the constant flare height control laws. All these differences were exolalned 
In detail In the orecedlng paragraphs and In Section 4. 


The block diagrams of Figures 4-4, 4-5 and the gain list of Table 4-II define 
the variable flare height system as tested In flioht on the Twin Otter. An 

V f 

unswwnetrical wind gain function fJO^lls shown In the upper left hand 

Vgo 

corner of Figure 4-4. This gain multiplies the closed looo sink rate and hence 
determines the closed looo flare time constant and -he slope of the closed looo 
line shown In Figure 4-fi. In the fast time simulation used for obtaining the 


evaluation results given in the following sections, •^hOF 


( CO ) 


* was Implemented 


Instead as a result of an error In the transmission of information from NASA to the 
contractor describing this control law. The asterisk denotes an upper limit of 


unitv for 


CO 


whereas the prime signifies a lower limit. Hence the gain with 


''GO 

the a arisk limits to a maximum of unity with headwinds rather than a minimum of 
unitv in tailwinds as Intended. The Implementation change reduces the closed looo 
time constant f si opes of the closed loop lines shown In Figure bv 35 percent 
for maximun 25 kt headwinds and 12 percent for maxlmim 10 kt tailwinds. The 
reduced slope results in a signal, to the control that the aircraft is too high at 
the flare entry. Note that the deterioration In landing performance will be 
greater for the spoiler cases, because for the spoiler control, there is no 
inhibiting time delay if the aircraft Is initially measured to he above the closed 
looo line as is the case for the pitch command fOR logic for FLR +1.5 sec in 
Figure 4-A) . 


5.2 0ETERHINI3TIC PERFORMANCE 

Simulation was used to obtain deterministic landing data for all six 
longitudinal control law configurations that were described 1n Section 5.1. Steady 
and shear winds vere Included as disturbances. Two tvpes of shearino winds were 
considered, one is linear and patterned according to FAA models (Reference 71 and 
the other is the sum of a logarithmic and a linear term. Both are described in 
Appendix A. The effect of headwinds uo to 25 knots and tailwinds uo to 10 knots 
were evaluated. (In the case of shear winds, these are the values at the reference 
height of 7.52 m, or 25 ft. above the ground whereas the wind magnitude at altitude 
is hioher). 
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Landinq time histories of the variable flare height control law with spoilers 
are given In Figures 5-lA and 5-lB. Three landings are shown - one with a standard 
FAA 25 knot shearing headwind, one with zero headwind and one with a standard FAA 
10 knot shearing tailwind. The last 30.5 m flOO ft) of each of the three 
approaches are shown. Rotation from the approach to the touchdown attitude Is 
performed through the flare. The app ach a^’tltude Is about -Q** for zero headwind 
(6® glide slope Is used), higher with the headwind and slightly lower with the 
tailwind. The touchdown attitude Is slightly positive for the headwind and zero 
wind cases and slightly negative with tailwind. (Touchdown attitude greater than 
-1® Is required for landing on the main gear first). Glide slope error Is less 
than 1 meter (3 ft) prior to the flare and during the flare the airplane deviates 
above the glide slope, as expected. Approach airspeed is 75 knots for the zero 
headwind case (A nominal approach speed of 71 knots Is used when the spoilers are 
not deployed). Pre-flare airspeed Is 73 knots for the headwind case and 76.5 knots 
for the tailwind as a result of the shears. The approach sink rate Is 4.15 m/sec 
(13.5 fps) with zero headwind, about 2.13 m/sec (7 fps) with the headwind and 4.73 
m/sec (15.5 fps) with the tailwind. Sink rate Is reduced through the flare to 0.76 
- 0.91 m/sec 2.5-3 fps) at touchdown. The sooller time histories Indicate a 
reduction of lift at the beginning of the flare, followed by an Increase of lift 
later on. The Initial lift reduction Is a result of the pitch flare height being 
higher than the commanded sink rate trajectory, (as sKB^'lf Figure 4-6) causing an 
Initially negative e', 45 seen in Figure 5-lB. The lift reduction with wind, 
particularly headwind, is more pronounced in the simulated system than in the 
system that was actually flown because of the difference that was explained in 
Section 5.1. The spoilers, throttles and elevator time histories here are given as 
deviations from their respective zero wind approach trim positions. The throttles 
are retarded to their flight Idle stops In all three cases shown. An uo elevator 
deflection of 12®-13® occurs through the flare to counter ground effect moments and 
provide the commanded rotation. Time histories of the major pitch command terms of 
the valrable flare height control law for the same three wind cases are given in 
Figure 5-lB. The block diagram of Figure 4-4 should be referred to In order to 
identifv the location of each of the variables shown. The Initial nose down 

command of 0^ does not propagate to because of the time delayed 
flare switch shown In Figure 4-4. The peak normal acceleration In the flare is 0.1 
q for the zero wind case and 0.16 g for the tailwind case. Additional landing time 
histories with the variable flare height control law are given In Appendix B. 

These time histories are for the no-sooHer configuration and for landings with log 
linear shears with or without spoilers. 











Figure 5- 1 A. Landing Time Histories; Variable Flare Height with Spoilers , 

5-12 0^^^-‘^',^2^V'- • 

(ft' 





Figure S>1B. Landing Time Histories; Variable Flare Height with Spoilen 
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Landing time histories of the low gain constant flare height configuration 
with spoilers are given in Figure 5-2 for the same wind shear cases as in Figure 
5-1. Rotation from the approach attitude to a slightly positive angle is accom- 
plished through the flare for all wind cases. Airspeed Is reduced through the 
flare from about 75 knots during the approach to 62 knots with headwind, 64 
knots with zero wind and 66 knots with tailwind. The targeted touchdown air- 
speed is 65 knots with the spoilers deployed (60 knots with the spoilers re- 
tracted). Sink rate is reduced from the approach value to 0.91 - 1.22 m/sec 
(3-4 fps) at touchdown. The spoilers reduce lift throughout the flare (be- 
cause the required rotation results in excessive lift). With zero wind, and with 
tailwind the throttles are being retarded all the way but In the headwind case 
the throttles are retarded only partially because airspeed is being actively 
controlled and more speed is being bled off In this case in comparison with the 
other wind cases, producing a speed arror that prevents the throttles from re- 
tarding all the way. Landing time histories of the other variants of the con- 
stant flare height system with the FAA and log-linear wind shears are given in 
Appendix B. 

Sink rate versus altitude trajectories are given in Figure 5-3 for the 
variable flare height control law without spoilers and in Figure 5-4 for the 
low gain constant flare height configuration without spoilers. As In the time 
histories, the last 30.5 m (100 ft) of the approach are shown. Each figure 
includes three trajectories, each of which was obtained with one of the follow- 
ing wind conditions: 25 knots log-linear headwind shear (as defined in Appendix 
A), zero headwind and 10 knots log linear tailwind shear. The variable flare 
height characteristic is evident in Figure 5-3 with sink rate reduction starting 
at about 15.2 m (50 ft) for the tailwind case and at 9.1 m (30 ft) for the head- 
wind case. With the constant flare height control law. Figure 5-4, rotation 
starts at 15.2 m (50 ft) for all wind cases but actual sink rate reduction starts 
at a slightly lower height for the tailwind case because of the higher initial 
sink rate involved. The touchdown sink rate dispersions are 0.67 - 1.43 m/sec 
(2. 2-4. 7. fps) for the variable flare height system and 0.66 - 1.26 m/sec (2.15 - 
4.15 fps) for the constant flare height control law. Additional flare sink-rate 
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trajectories, for the various control laws and deterministic wind conditions that 
were evaluated, are given in Aopendix B. 

The variation of touchdown sink rate, range and pitch attitude with wind 
magnitude is shown in Figure 5-5 through 5-10. Each one of the six figures is for 
one of the six control law configurations that have been evaluated. There are 
three curves in each olot - one for steady wind, one for the linear FAA shears and 
one for the log- linear shears (all wind models are defined in Appendix A). Wind 
magnitudes vary from 25 knots headwind to 10 knots tailwind. No stochastic 
disturbances, i.e., turbulence, were included. Deterministic touchdown results for 
the variable flare height control law without spoilers are given in Figure 5-5. 

The log-linear shears cause the largest dispersions on all three touchdown 
variables. This is due to the fact that with this tvpe of shear the variation of 
wind close to the ground is more rapid than w*th the FAA linear shears. The 
hardest touchdown sink rate is 1.^0 m/sec (4.5 fps) with 25 knots of loo linear 
headwind shear and the softest is 0.57 m/sec (2.2 fosl with 10 knots of log linear 
tailwind shear. Sink rate variation with the FAA shears and steady winds is very 
small. Touchdown range is measured down the runway referenced to the Glide Path 
Intercept Point (GPIPl. The shortest touchdown distance for the variable flare 
height configuration without spoilers is 6.1 m (20 ft) beyond the GRIP with 25 
knots log linear headwind shear, as seen In Figure 5-5. The 10 knots log linear 
tailwind shear produces the longest landing for this configuration at 88.4 m (290 
ft) beyond the GRIP. Pitch attitude varies between -0.4® and 1.9®. The equivalent 
results for the variable flare height configuration with spoilers are given in 
Figure 5-6. Sink rate dispersion Is about the same as without the spoilers, range 
dispe'^sion Is slightly improved and touchdown attitude with spoilers varies from 
-0.7® to 1.4®. The tendency towards lower attitude with spoilers seems to be 
associated with the spoiler's tendency to Increase lift at the end of the flare 
(see Figure 5-lA) and it could also be related to higher airspeed used with the 
spoilers. Deterministic touchdown results for the low oain constant flare height 
control law without spoilers are given In Figure 5-7 and with spoilers In Figure 
5-8. Here again, the log-llnsar shears result In the largest dispersions. Sink 
rate with no wind is significantly harder, at 1.22 m/sec (4.0 fos), for the 
constant flare height configuration in comoarlson to 0.85 m/sec (2.8 fps) for the 
variable flare height control law. (Note that a harder touchdown sink rate will 
tend to reduce touchdown dispersions). Sink rate with the constant flare height 
does not Increase significantly with headwinds. Sink ■ ..^.e dispersion is slightly 
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better with the constant flare height. Range has the same oeneral tendency of 
being short with headwind and long with tailwind. Range dispersion Is slightly 
Improved with the constant flare height configuration In comparison to the variable 
flare height. The variation of touchdown pitch attitude with constant flare height 
Is much smaller than with the variable flare height. The use of spoilers with the 
low gains constant flare height configuration results In an almost Insignificant 
Improvement of sink rate control, but range dispersion Is reduced from 80.8 m (265 
ft) without spoilers to 57.9 m (190 ft) with spoilers. The use of spoilers also 
result, with this configuration. In a slightly higher touchdown attitude. This Is 

because with this configuration the spoilers reduce lift throughout the flare. The 

determini si tic touchdown results for the high gain constant flare height 
configuration, without and with spoilers, are given In Figures 5»9 and 5<>10, 
respectively. Sink rate control with high gains and no spoilers Is Improved In 
comparison to the low gain configuration but with spoilers It Is about the same. 
Range dispersion Is reduced for the high gain configuration. It Is 65.5 m (215 ft) 

without spoilers and 27.4 m (90 ft) with spoilers. Touchdown pitch attitude with 

the high gains Is similar to that obtained with the low gains. 

Table 5-V sunmarlzes the results given In Figures 5-5 through 5-10 (Touchdown 
airspeed Is an additional variable Included In the table). The numbers given In 
the table are the extreme values of each variable for all wind conditions that were 
evaluated. The hardest sink rate values are 1.40 m/sec (4.6 fps) with the variable 
flare height control law and 1.25 m/sec (4.1 fps) with the constant flare height 
configurations. The low gain constant flare height conf louratlon with spoilers and 
the high gain constant flare height configuration without spoilers provide less 
sink rate dispersion than the other configurations. The use of spoilers has no 
significant affect on sink rate control with the variable flare height control 
law. The use of spolle-s reduces sink rate dispersion for the low gain constant 
flare height control law but not for the high gain version. Touchdown range 
dispersion is not Improved by any significant amount by the use of spoilers with 
the variable flare height control law. Range control of tne low gains constant 
flare height system without spoilers 1 ; about the same as with the variable flare 
height control law. With spoilers, however, the low gain constant flare height 
system obtains a 30 percent Improvement in range control for these deterministic 
disturbances. Range control with the high gain constant flare height system Is 
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Figure 5-7. Effect of Deterministic Wind on Touchuown,Constant Flare Height, 
Low Cains, No Spoilers 
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Figure 5-9. Effect of Detenninistic Wind on Touchdown, Constant Flare Height, 
High Gains With Spoilers 
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TABLE 5-VA DETERMINISTIC TOUCHDOWN NUMMARY 
(METRIC UNITS) 



DISTURBANCES ARE STEADY WINDS. FAA SHEARS AND LOG LINEAR SHEARS VARYING FROM 


TABLE 5-VB DETERMINISTIC TOUCHDOWN SUMMARY 
{ENGLISH UNITS) 



FROM 25 kt HW TO 10 kt TW. 


improved by 20 percent in comparison to the low gain version. The use of spoilers 
results in a further improvement of range control by 60 percent. Thus, the high 
gain constant flare height system with spoilers produces far better range control 
than the other configurations. The touchdown attitude results of Tabe 5-V indicate 
a somewhat uncomfortably low attitude for all configurations as an attitude greater 
than -1* is required to assure landing n the main gears first. The variable flare 
height with spoilers is worst in this respect at -0.7* and the constant flare 
height configurations with spoilers are the best, touching down at a minimum 
attitude of 0.2* or 0.3* , The targeted touchdown airspeed is 65 knots with the 
spoilers deployed and 60 knots with the spoilers retracted. The higher airspeed is 
used with spoilers to comoensate for an increase in stall speed. The minimum 
touchdown speed for all spoiler configurations is 6.0 to 6.5 knots slow and for the 
no-spoiler configurations is 4.0 to 5.5 knots slow. The constant flare height 
configurations touchdown at about 1.5 knots slower than the variable flare height. 
All the minimum touchdown airspeeds are associated with the 25 knot log-linear 
headwind shear whereas other wind conditions result in higher touchdown airspeeds. 

To summarize the results with deterministic winds, the low gain constant 
flare height control law provides somewhat better sink rate, range and attitude 
control than the variable flare height configuration. The high gain constant flare 
height system with spoilers provides much better range control than all the other 
configurations. The. log-linear shear is the most demanding of the winds evaluated. 

5.3 STOCHASTIC PERFORMANCE 


The simulation that had been developed for this study was used to evaluate 
landing performance of the six longitudinal control law configurations with 
stochastic as well as deterministic disturbances. Horizontal turbulence, vertical 
turbulence and MLS beam noise according to the models that are defined in Appendix 
A were included in the simulation. Simulation runs were made with limiting FAA 
shear winds, moderate turbulence and beam noise. Probability distributions of the 
major touchdown variables were obtained and the results are summarized in this 
section. Data with respect to activity of the controlled variables and controllers 
on the glide slope with turbulence and beam noise were also obtained and are 
summarized here. Landing performance was also evaluated as a function of wind and 
turbulence level and the results are sumnarized in this section. Note that the 
limiting wind and turbulence levels are higher than the average levels that would 
be encountered in actual flight. While results using these levels provide a good 


measure of the effectiveness of the control system heinq evaluated, thev are a 
conservative indication of the abilitv to meet statistical criteria which are based 
on smaller average levels of disturbances. 

5.3.1 LANDING PERFORMANCE WITH LIMITING WINDS 

Table 5-VI is a summary of the longitudinal landing performance with limiting 
shearing winds, moderate turbulence and beam noise. The left most column of this 
table defines the oerformance goal that was used in this study. As explained in 
Reference 1, requirements for automatic landing systems of CTOL transport aircraft 
have been established but not for the STOL airplane. Consequently, performance 
goals had to be defined for this study, based on CTOL requirements scaled down to 
the STOL runway and landing geometry and taking into account geometrical and 
physical characteristics of the Twin Otter airplane. The mean taraeted sink rate 
was selected to produce a comfortably soft landing. The two sigma land hard is a 
design objective which would produce acceptable sink rate control and a 10“^ 
hazard probability is attached to the exceedance of the Twin Otter airplane's gear 
strength. On range, the goal for the mean was computed assuming the airplane to he 
on the glide slope at flare initiation and executing the commanded h/^ flare 
trajectory with no deviations. Two sigma dispersions of m f ±200 ft) were 
scaled down from CTOL requirements. The 10”^ land short reauirement provides for 
landing within the STOL runway's safety underrun area f Figure 3-1). The 10-^ 
land long requirement depends on the airplane's stopping distance and the runwav 
length. Such considerations were outside the scope of this study and therefore the 
mmber shown in Table 5-VI is simply a linear extrapolation of the two-sigma 
dispersion to the 10*^ probability level. The 10-^ low attitude goal is based 
on the hazard of hitting the ground with the nose wheel first. All the control 
laws that were evaluated tended to produce flat touchdowns and therefore upper 
limits on attitude were not required. 

The actual performance results that are given in Table 5-VI were comouted 
assuming a 70 percent probability of encountering a 25 knot shearing headwind and 
30 percent orobabilitv for a 10 knot shearing tailwind. The 70/30 split is based 
on the results of a survey of 79 major U.S. airports looking at runwav orientation 
wit^" respect to the prevailing winds. This ratio is apoli cable for airports that 
have bi-directional landing aids on the maii'^r runwav. The assumption that each 
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TABLE 5-VIA LONGITUDINAL PERFORMANCE SU»«ARY (METRIC UNITS) 

( 70 % HW, 30% TW) 







D I STURBANCES * 

70% PROBABILITY FOR 25 kt SHEARING HW, 3.73 kt 
30% PROBABILITY FOR 10 kt SHEARING TW. o = 2.67 kt 


TABLE 5-VI B LONGITUDINAL PCrtFORMANCE SUMMARY (ENGLISH UNITS) 

(701 HW, 30* TW) 



landing Is made with either a limiting headwind or a limiting tailwind Is 
conservative since the probability of encountering 25 knot headwinds or 10 knot 
tailwinds Is significantly lower according to Reference 7. Simulation runs were 
made with 25 knot shearing headwind and 3.7 knot RMS horizontal turbulence or with 
10 knots shearing tailwind and 2.7 knots RMS horizontal turbulence. All runs were 
made with 1.5 knots RMS vertical turbulence. Shears according to the FAA models 
only were used In obtaining the stochastic data. 

The results of Table 5-VI Indicate that the 10"® land hard requirement Is 
exceeded by the variable flare height control law with or without spoilers. All 
constant flare height configurations provide acceptable to very good sink rate 
control. The low gain system without spoilers Is the poorest of the four in this 
respect ^KJt Its performance Is very significantly Improved with the use of the OLC 
spoilers. Sink rate control with the high ofin constant flare height control law 
Is good even without the use of spoilers and It Improves somewhat with spoilers. 

All control law configurations meet the 10"^ land short requirement but the 
high gain configuration with spoilers is the only one to meet the 10"*^ land long 
requirement. For the variable flare height control law, the use of spoilers 
improves range control by about 20 percent. Range control with the low gain 
constant flare height configuration without spoilers is about eguivalent to that of 
the variable flare height without spoilers but the use of spoilers here falls to 
produce a significant improvement. It is suspected that switching of gains when 
transitioning from glide slope track to flare, adversely affects the range control 
with the low gain constant flare height configuration and therefore the improved 
sink rate control with sooilers does not translate to an eguivalent improvement in 
range control. Range control with the high gain constant flare height system 
without spoilers is as good as with the variable flare height control law with 
sooilers. Excellent range control Is obtained with the high gain constant flare 
height configuration with spoilers. The use of spoilers reduces range dispersion 
by 40 percent for this configuration. 

The high gain constant flare height configuration with sooilers is the only 
one that meets the 10”6 low attitude requirement. All the configurations that 
were evaluated produced rather flat touchdown attitudes. Attempting to arbiti^irilv 
command a higher touchdown attitude with this light wing loading airplane, would 
result in a tendency towards long landings. 
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The nose down tendency was aggravated by the 4 knot Increase in approach 
soeed to account for a higher stall speed with spoilers deployed to their nominal 
value. If an automatic system were installed to close the spoilers for a downward 
gust O'* shear correction, the spoilers closed stall soeed may be approoriate as the 
basis for establishing approach airspeed. When spoilers are used a lower approach 
speed would increase the approach pitch attitude and probably reduce the touchdown 
dispersions. 

Deviations of more than 3.66 m (12 ft) from the glide slope at a decision 
height (Ah window) of 30.5 m (100 ft) would result in a reouirement to abort the 
approach. The automatic landing system is reguired to be within this window on a 
two sigma basis. The variable flare height configuration without soo11«*rs deviates 
somewhat from this requirement on the low side. (The asymmetrical nature of the 
deviations from the glide slope result from the dominant effect of the shearing 
headwind). All other configurations meet this glide slope window requirement. The 
low gain constant flare height system tracks the glide slope better than the 
yariable flare height control law and the high gain version produces better glide 
slope tracking performance than the low gain version. The use of spoilers improves 
glide slope tracking performance with all configurations. 

The probability distribution plots from which the data in Table 5-VI was 
extracted, are given in Appendix B. 

5.3.2 CONTROL AaiVITY 


Activity data for controlled variables and controllers on the glide slope 
with moderate turbulence and MLS beam noise are summarized in Tables 5-VII and 
5-VIII. Data were obtained from simulation and are given for all six evaluated 
control law configuraions. The results of Table 5-VII are with moderate turbulence 
and MLS beatn noise, as defined in Appendix A. The results of Table 5-VIII are with 
turbulence only. 

Table 5-VII indicates that glide slope track accuracy with the low gain 
constant flare heigiit control law is s’ightly better than with the variable flare 
height configuration. The high gain configuration provides a glide slope tracking 
accuracy that is significantly better than with the other configuration. Glide 
slope tracking accuracy is improved by the use of the OLC spoilers. The 
improvement Is significant for the variable flare height and the low gain svstem 
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3.73 kt, 0 » 1.5 kt, » 0.06“ 0 « 152 m (500 ft) 



3.73 kt 


but fl>1nor for the hlqh qain version. Sink rate and normal acceleration activity 
are about equivalent with a11 conf 1 quratl ons but the low qain control law has 
sllqhtly hlqher sink rate activity and the hlqh qain system has sllqhtly lower sink 
rate activity than the other cnnfl quratl ons. Alrsoeed control Is also about 
equivalent for all conf 1 quratl nns with the high qain configuration being sllqhtlv 
better than the others. Pitch attitude and elevator activity with the constant 
flare height control laws are somewhat hlqher than with the variable flare height 
system. Throttle activity with the constant flare height configurations Is almost 
twice as high as with the variable flare height. This Is a result of not using 
longitudinal acceleration feedback to the throttle in the constant flare height 
systems (see Table 5-III). Glide slope spoiler activity Is 8® 'IMS with the 
variables flare height, somewhat higher (10®) with the low gain constant flare 
height system and much lower (4.8®) with the high gain configuration uS a result of 
the higher pitch gains. The use of spoilers Improves glide slope tracking accuracy 
and reduces sink rate, normal acceleration, airspeed, pitch attitude, elevator and 
throttle activities. 

Table 5-VIII sunmarizes glide slooe activity with horizontal and vertical 
turbulence but no beam noise. The difference between the results of Tables S-VII 
and 5-VIII Is a result of US beam noise being Included In the former. Glide slope 
tracking accurajfy^ithout beam noise looks much better because Ah^ is raw glide 
slope deviation filtered with a 0.3 second first order lag, such that beam noise 
shows up directly on this variable. Without the beam noise, the improved glide 
slope tracking accuracy and sink -rate control with the high qain configuration are 
more apparent. Beam noise contributes some sink rate activitv hut little normal 
acceleration. It has a small impact on alrsoeed, pitch, elevator and throttle 
activities. Beam noise contributes to spoiler activity only with the low gain 
constant flare height control law. 

The glide slooe tracking accuracy ^iven in Table 5-VIII is a good indication 
of anticipated accuracy at the 30.5 m (100 ft) decision height because the 
contribution of beam noise at 30.5 m (100 ft) is much smaller than at 152.4 m (50P 
ft) as a result of the proximity to the transmitting autenna. The PMS results of 
Table 5-VIII can be doubled to obtain an estimate of the two s^gma excursions. 

Glide slope tracking accuracy is thus ±3.90 m f±i?.8 ft), on a two sigfla basis, for 
the variable flare height and +3.66 m (12.0 ft) for the low qain constant flare 
configuration, both without spoilers. This is margina’lv acceptable comparison 
to the required accuracy of t3.66 m (±12 ft). With spoilers, glide slooe trackino 
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accuracy of these two configurations is good. The performance of the high gain is 
excellent with, or without, spoilers. 


5.3.3 LANDING PERFORMANCE WITH REDUCED WIND LEVELS I 

i 

I 

Landing perform,* nee results with limiting winds and moderate turbulence were 
given ‘n Section 5.3.1. Results with reduced wind and turbulence levels are given 
and discussed here. Mean and two sign? touchdown results as a function of wind and 
turbulence level are given in Figure 5-11 for the variable flare height 
configuration without spoilers and in Fipure 5-12 with spoilers. Data were taken 
with FAA wind shears of 25, 12.5, and 0 knots headwind and 5 and 10 knots 
tailwind. The magnitude of the horizontal turbulence w?s varied with the 
deterministic wind level as shown in the figures. A constant level of 1.5 knots 
vertical turbulsnce and beam noise were included in all cases. No ho»*izontal 
turbulence was included with the zero headwind cases and therefore dispersions at 
this point are due to vertical turbulence and beam noise only, and vertical 
turbulence is the major contributor. 

The resuHs of Figure 5-11 indicate that the mean touchdown sink rate is 
essentially constant regardless of the wind magnitude. The two sigma land hard 
value is about 0.30 m/sec (1.0 fps) higher than vhe mean sink ra-W!*with zero wind 
(and horizontal turbulence). The worst case for sink rate control is at 25 knots 
headwind and the associated 3.73 knots of horizontal turbulence. The two signa 
land hard value there is 1.83 m/sec (5.0 fps), or 0.88 m/sec (3.2 fps) harder than 
th' mean sink rate. The mean touchdown range varies from 55.4 m fi85 ft) with 25 
knots shearing headwind to 88.4 m (280 ft) with 10 knots shearing tailwioj. fRanae 
is measured vnth respect to the GRIP). Touchdown range dispersion, from two si am a 
chort to two signa long, varies from 45.7 m (15C ft) at zero wind to 158.5 m (520 
ft) with headwind and 3.73 knots of horiz:;ntal turbulence. The dispersion is 
smaller with tailwind (because of the smaller amount of turbulencel . The mean 
touchdown pitch at*itude is fairly constant at about i'* for all headwind magniirj'ies 
but it drops to ' • • ..ith the 10 knots shearing tailwind. Most of the attitu;'£ 
dispersion is caused by the vertical turbulence. It is 1.2° from the mean to the 
two cigno low attitude with zero headwird and 1.55° for the maximum headwind or 
tailwind. The lowest attitude on a two sigi^a basis is -1.55° with the 10 knots 
tailwind. Glide slope tracking dispersion at the 30.5 m (100 ftl decision height 
varies from 1,83 m .6 ft) with vertica) turbulence and beam noise onlv, to 7.6? m 
(25 ft) with the full 25 knots headwind and the associated turbulence. The results 
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shew a -0.76 m (-2.5 ft) mean nevlation from the glide slope with zero wind. This 
is a result of an offset in the simulation affecting the beam complementary filter 
and resulting in a small constant beam error even in the absence of disturbences. 
The mean tends to be lower with the shearing headwinds and higher with the 
tailwinds. 

The results for the same configuration with spoilers are given in Figure 
5-12. All the dispersions are significantly reduced. Sink rate control is 
improved by more than 40 percent, range by 60 percent, attitude by more than 50 
percent and glide slope track by 30 percent. The variation of the mean touchdown 
range and attitude as a function of wind has increased, however. 

To sutmarize, all dispersions grow significantly with the magnitude of 
horizontal turbulence. Pitch attitude is the only touchdown variable evaluated 
that was more affected by vertical than by horizontal turbulance. The use of 
spoilers with this variable flare height control law produces a very marked 
improvement in the control of range, sink rate and glide slope deviations in the 
two sigma region. 

The probability distribution curves, on which the data in Figures 5-11 and 
5-12 are based, are given in Aopenoix B. 
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6.0 LOCALIZER TRACK AND RUNWY ALI6Nf€NT 

For good lateral automatic landing prformance, the runway alignment maneuver 
must be Initiated with the aircraft stabilized on beam certerllne. Thus, the 
design of a stable, tight localizer track mode is a necessary orerequislte to 
defining an advanced automatic decrab control system. The assumotlon is made that 
even close in curved path approaches are selected so that the capture algorithm 
will place the aircraft on the localizer sufficiently early to allow transients to 
settle. Since the capture will have minimal Impact on landing performance. It Is 
not considered in this study. For Category IIIA landings, automatic roll-out 
guidance is not required, and it also is not included in this study. 

During the design of the lateral landing system for the Augmentor Wing STOL 
vehicle. Reference 3, extensive tradeoff and optimization studies were performed to 
define the recommended forward slip algorithms. With few exceptions, the Twin 
Otter lateral directional dynamics and system requirements are very similar to the 
Augmentor Wing. Thus, maximun useage was made of the tradeoffs conducted during 
the previous work to provide a sound basis for this Advanced Autoland lateral 
landing study. This allowed a more extensive evaluation of the Twin Otter 
peculiarities, with tradeoffs and optimization directed specifically toward the 
problems peculiar to a light wing loading STOL aircraft. 

This section describes the performance of the recommended localizer track and 
runway alignment control laws that are described in Section 4. Performance was 
evaluated through the use of simulation with deterministic and stochastic 
disturbances. The effects of sensor errors, trim changes and svstem variations 
were evaluated and the results are given and discussed here. Control law 
variations that were studied in the process of defining the recommended control law 
are also discussed in this section, as well as failure effects and system 
limitations. 

6.1 PERFORMANCE 

A performance summary of the recommended lateral landing system as described 
in Section 4 is presented here, with the detailed supporting data included in 
Appendix C. 

The following subjects are presented and discussed here: 

1. Landing time histories. 

2 . Statistical landing performance in stochastic disturbances only. 

3. Landing performance over the total landing environment. 

4. Deterministic variations and off-nominal conditions. 

5. Localizer track activity in stochastic disturbances. 
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Landing Time Histories 

Landing time histories, as obtained from simulation, with the recommended 
control law are given In Figure 6-1 (A and B). Six time histories are shown, 
starting at a C.G. height of 91.4 m (300 ft) above the ground and terminating at 
touchdown. The flrs^ three are all with 15 knots right shearing crosswind and the 
first Is with 10 knots shearing tailwind, the second Is with zero headwind and the 
third Is with 25 knots shearing headwind. The second group of three landings Is 
with 15 knots right steady crosswind and the same headwinds. Wind shears are 
patterned after FAA models as described in Appendix A. The Initial crab angle Is 
20° with the shearing crosswind and IZ"* with the steady crosswind because the 
crosswind magnitude is higher at altitude with the shear. At touchdown, the 
airplane is aligned with the runway heading , within one degree, for all wind 
cases. There are no heading overshoots. Crosstrack velocity at touchdown is less 
than -O.llA m/sec (-0.375 fps away from the windl. hank angle at touchdown Is less 
than 5° (Into the wind) for all cases. The control law does not have a wing 
leveling mode. The sideslip Induced lateral acceleration is about -0.087 a for all 
wind cases. Ten degrees of rudder are needed to maintain runway heading and 
control wheel excursions of less than 10° are used. The crosstrack displacement at 
touchdown varies from 0.61 m (2 ft) to the right to 1.07 m f3.5 ft) to the left. 
Additional variables for the same six landings are shown In Figure 6-iB. The ”CG 
traces clearly show the Impact of the headwind variations. The alignment model, 

’l>H|, as defined in Figure 4-7 starts out being equal to the crab angle above 45.7 
m (150 ft) and It Is reduced to zero at 15.2 m (50 ft). The rudder alignment 
commetfid, 5C alI®. (which Is the .rudder command during alignment, 
excluding ?he yaw rate path) hits the 16° limit x Kn_, Figure 4-7) for 

shearing crosswind cases. The 6^ trace Is the total rudder command, 
including the yaw rate path. Rudder command does not exceed 11° for any landing. 

S tatistical Landing Performance 

The landing performance of the recommended lateral landing system Is 
surmarized in Table 5-i both for stochastic disturbances only and over the total 
landing environment including the deterministic disturbances defined in Table 
6-1 1. These data were obtained with limiting atmospheric and MLS disturbance 
levels, with average limiting longitudinal winds and shears, without accounting for 
the low occurrence probability of tnese disturbance levels, and thus it represnets 
a conservative estimate of landing performance. Moreover, Crows Landing Stoloort 
(the airport at which the Twin Otter was tested! geometry, with localizer antenna 
placement of 1356 m (a450 ft) with respect to the elevation antenna location was 



t 
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used to give la ■ 2.91 m (9.55 ft) lateral beam noise level at an altitude of 
30.5 m (100 ft). In an actual 5T0L port Installation a number less than 1355 m 
(4450 ft) can probably be established. 

The lateral performance goals are alven In the first column of Table 5-1. 
Performance requirements for Cateoorv III automatic landing systems have been 
defined for CTOL transports but not for STOL aircraft, as discussed In Reference 
1. Thus, design goals had to be defined for this study. Most of these goals were 
specified by the NASA Statement of Work (SOW) or are determined bv vehicle 
constraints. The allowed lateral two sigma landing dispersion was specified as 
±4.57 m (15 ft) In the SOW. However, the FAA AC20-57A (Reference 7) requirement 
for not landing any closer than five feet to the runway edge, effectively specifies 
that for the Twin Otter vehicle on a 30.5 m (100 ft) runway, the probabilltv of 
exceeding 11.88 m (39.0 ft) of lateral deviation at touchdown must be less than 
10"®, which Is equivalent to a two sigma dispersion of ±4.86 m (15.9 ft). In 
Table 5-1, this 10“® value Is being used along with the SOW requirement for two 
sigma lateral dispersion. Applying the above 11.88 m (39.0 ftl constraint to 
lateral deviation during rollout, the touchdown heading error becomes 10.8®. This 
Is based on the assumption that lateral rate just after touchdown is approximately 

equal to 11 ^/ 57.3 times the touchdown course error, and that a correcting lateral 
acceleration of 0.2 g's Is applied to arrest the lateral velocity. The 

corresDondinq 2 sigma value for ^5 4.40 degrees. The touchdown bank angle 
requirement is based on aircraft geometry. Wlno contact for the Twin Otter 
airplane occu-s at 19® of bank. The touchdown crosstrack velocity limit is 
normally dependent on touchdown bank anale. For conservatism, a one oea*" touchdown 
lateral velocity limit of 3.048 m/s (10 fos) was considered as the 10-® 
requi'*ement. The results of Table 5-1 indicate that the landing oerformance 
satisfied all the design goals by wide marains. The probability distributions for 
the landing parameters of Table 5-1 are given in Appendix C. 

Deterministic Variations 

Many deterministic variations were considered during this study. Those which 
significantly impact lateral landing performance are listed in Table 5-II along 
with their effect on the aircraft performance Indication parameters, with a more 
complete discussion given in Lhe next section. 


TABLE 6-1 LATERAL PERFORMANCE SUMMARY 


PABAMETER 

GOAL 

ACTUAL 

i 



STOCHASTIC ONLY 

TOTAL POPULATION 

TO 

U • 

-0.365 (-1.2) 

-0.365 

(-1.2) 

m (ft) 

20 <4.57 (15.0) 

-3.35 (-11.0) 

-4.419 

(-14.5) 


10‘® <11.88 (39.0) 

-•^^.58 (-21.6) 

-8.9 

(-29.2) 

TO 

u - 

-0.122 (-0.4) 

-0.122 

(-0.4) 

m/sec (fps) 

2o - 

±0.518 (±1.7) 

±0.518 

(±1.7) 1 


10“® <3.048 (10.0) 

-1.37 (-4.5) 

-1.37 

(-4.5) j 


u - 1 

0.0 

0.0 

1 

i 


2c . 

' r+2.2 ! 

r +2.4 

i 

} 


1 (_-0.8 i 

i-..3 

» 

\ 


lO"® <10.8 1 

+7.2 i 

i 

+7.7 

i 

*TD • 

u - 

+4.7 i 

1 

+4.7 


1 i 

20 - i 

• \ 

±2.0 : 

to 

• 

i 

! i 

1 

-6 



' 

! 

1 

j 

10 <19.0 1 

-I-9.4 

+9.8 


^WINDOW ; 

i 

u i 

1 

+0.457 (+1.5) 

+0.457 

(+1.5) 

m (ft) 

1 

1 

2o (±25.0) 1 

±3.81 (±12.5) 

±4,11 

(±13.4) : 


NOTES: 


1. This table defines both stochastic only and total population landing 
performance, where the latter includes deterministic disturbance effects. 

2. Only those requirements specified in the SOW or dictated by aircraft/ 
geometr^^ limitations are listed. 

3. For the lateral axis, means should all be zero when evaluated over the 
total environment. Since actual results are presented for liiLit. rg wind 
from thv right, significant mean values are obtained. 

4. is lateral cracking error at the 30.48m (100 ft) approach window. 

WINDOW 

5. For conservatism, results are presented for limiting crosjwind, shear, tur- 
bulence, and MLS noise levels for averaged limiting headwind and tailwind 
conditions, without accounting for the "-^4% probability occurrence of 
these limiting atmospheric disturbances. 

6. Deterministic variations included in total population aistributions are 
MLS azimuch bias, course datum error, accelerometer gyro errors, aiid 
approach speed variations. 
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TABLE 6- 1 1 LATERAL DETERMINISTIC DISTURBANCE EFFECTS 



These results are graphically presented in Appendix C. 

See discussion of Rudder Trim Shift in paragraphs 6.2.3 and 6.3.2. 


The performance Impact of these deterministic variations is relatively 
small compared to atmospheric disturbance induced landing dispersions. However, 
the lateral displacement error at the intended touchdown point due to the 0.0935*’ 
one sigma azimuth bias is significant. 

Activity in Disturbances 

A summary of aircraft state variations and control activity as obtained from 
simulation with limiting turbulence levels and beam noise is given in Table 6-III. 
All activity levels are acceptable. 


TABLE 6-III LOCALIZER RMS ACTIVITY 

" — % I— ' . - ■ 


♦ 

deg 


0.85 

♦ 

deg 


1.60 

y 

m/s (fps) 

0.427 

(1.40) 

*y 

9 


0.010 

• 

♦ 

deg/s 


0.68 

«R 

deg 


0.80 

5 

w 

deg 


2.0 

y 

m (ft) 

1.34 

(4.40) 

NOTE 

• ‘^BN 

» 2.9 m (9.5 

<^sG - 


In summary, the recommend''d system provides landing performance compatible 
with Cat IIIA landing requirements. 

6.2 SENSOR ERRORS AND SYSTEM VARIATIONS 


The following sensor and system errors and off-nominal conditions were 
evaluated as part of this study and are discussed here: 

1. MLS guidance system errors 

2. Sensor errors 

3. Rudder trim changes 
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4. Sensitivity analysis 

5. MOOILS discretization effects 

6.2.1 MLS ERRORS 


For a Cat IIIA lateral automatic landing system, o'lly OME and azl.nuth errors 
can Impact landing performance. The effects of OME bias Is small, and only azimuth 
variations need be considered. Assuming that the azimuth antenna Is located 6C9.6 
m (2000 ft) from the nominal touchdown point, the 0.0935* one sigma bias level 
directly yields ±0.99 m (±3.26 ft) lo lateral touchdown variations; and 4.88 m 
(16.0 ft) If extrapolated linearly to the 10*6 (4.9o) probability level. Thus, 
this term {.lone Is as significant as limiting wind levels on lateral touchdown 
dispersion. To limit the Impact of this bias In the low probabllltv region, a 0.3* 
3o cutoff level was used, since It appears reasonable that the near field beam 
monitor threshold should be no greater than 0.3*. With this assumption, the 
lateral deviation due to stochestic and deterministic disturbances at the 10"® 
probability level is 8.9 m (29.2 ft), as shown in Table 6-1. 

6.2.2 SENSOR ERRORS 


The three sensor on board the aircraft which contribute significantly to 
landing performance are the compass system, the lateral accelerometer, and the 
vertical gyro. 

Since the compass system signal provides the main input to the align 
maneuver, its errors are significant. The expected accuracy is ±4.0 degrees on a 
4.5o basis. As expected. It vields touchdown misalignment on a one to one basis, 
with relatively small Impact on other landing parmneters. 

The runway axis lateral acceleration signal Is used both for navigation 
filter au^entatlon and wing down compensation, and its error characteristics are 
of a dynamic nature. An equivalent offset of ±0.014 g, 4.5o (or 0.0031 q, lo 
as shown in Table 6-II) was assumed. The induced landing variations are very 
dependent on landing systc-n design; for the recomnended configuration with washed 
out acceleration, this acceleration inaccuracy yields 0.45 m (1.48 ft) one sigma 
lateral dispersion as discissed in detail in Apaendix C. 
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The roll attitude vertical gyro errors are also dynamic In nature, a 
combi natl »on of steady errors and acceleration Induced erection and drift errors as 
discussed In Appendix C. The *0.35® equivalent lo error has about the same 
landing performance Impact as the accelerometer (see Table 6-II). 

6.2.3 RUDDER TRIM (31ANGES 


The usage of the rudder command signal for crossfeed ourrioses, has an Impact 
on the landing accuracy as discussed In paragraph 6,3.2. The effect Is dynamic ^n 
nature and Is a function of the glide slope capture maneuver characteristics and 
the altitude at which this maneuver 1s executed. A ±1.5® rudder deflection on a 
3<j basis was determined fas showr, In oaragraoh 6,3.21 as the effective erro»* 
Introduced via this oath. The impact on lateral touchdown deviation Is ±1.07 m 
(±3.5 ftl three slc^a. 

6.2.4 SENSITIVITY ANALYSIS 


A sensitivitv analysis study was conducted to examine thf effect of gain and 
time constant variation on system oerformance. The analysis Included t’te variation 
of most of the gains bv ±6 db around nominal, with the exceotlon of the yaw to roll 
crossfeed gains for which the variation was ±20 percent. The landing 
histories are given in Appendix C. As expected, none of the variables has a 

pronounced effect, with and i 2 having similar impact as crossfeed 
parameters variations. 

6,2.5 MDDILS DISCRETIZATION EFFEaS 


MODILS is an experimental microwave landing aid that used initially in 
flight tests. It has a narrow azimuth coverage and a qra;;uia- u'imuth signal. The 
impact of MODILS finite resolution (refer to Appendix Al was investigated to 
establish the feasibility of using higher guidance proportional loop gain (Kyi. 

The MODILS beam resolution of 0.1®, renders 4.72 m fl5.5 ftl lateral deviation 
discretization effect at an -.iltitude of 274 m (900 ftl. Localizer track 
oerformance with tlie discretization Incorporated into the simulation is oiven as 
Kiaure 6-2 where a rudder ramo going from zero to 1.5® in 10 seconds was used as 
the disturbance. The "’suits art optimistic since a discretization level of ±3.28 
m (±10 ftl was used. Note that for a doubled proportional gain, absolulelv no 
lateral deviation performance improvement is gained. On the contrarv, a higher 
limit rvcle freouencv is observed. 
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The resolution of the HLS guidance systen Is 0.01* In azimuth, an order of 
magnitude better than MOOILS, yielding 0.457 m (1.5 ft) steps at an altitude of 274 
m (900 ft). The scaling of the digital autopilot Is 1.22 m/bIt (4 ft/bit) with 
MOOILS but It has been changed to 0.0762 m/bIt (0.25 ft/bit) for the final approach 
phase with MLS guidance. Thus, the MLS beam Information can be considered 
continuous. 

Figure 6-3 sunmarlzes the effect of MOOILS resolution on statistical landing 
performance. Results are given for the discretized beam signal and for a 
continuous signal. Data were taken with turbulence and beam noise and repeated 
with turbulence only because beam noise acts as a dither, masking the effects of 
the poor resolution. This is most obvious on y^p^ ^he beam noise the 

results for the discrete and the continuous signal are very similar. Without the 
noise, the continuous beam yields a much smaller lateral dispersion than the 
discretized signal. The sane effect Is apparent to a smaller degree on the 
touchdown bank angle. 

6.3 SYSTEM OPTIMIZATIOM 


This section describes the results of tradeoff studies that were carried out 
in order to arrive at the recommended localizer track and runway alignment control 
laws that are described In Section 4. Extensive tradeoff studies were made In the 
course of a previous automatic landing study for the Augmentor Wing airplane which 
is a powered lift experimental STOL airplane. Some of the major results of that 
study are considered to be directly applicable to the Twin Otter. They are listed 
here and desc»‘1bed in detail in Reference 3. This section concentracts on 
presenting the results of tradeoff studies that were directed at problems peculiar 
to the Twin Otter as a representative light wing loadinq STOL airplane. 

6.3.1 RESULTS BASED CM THE AUGMENTOR WING STUDY 


The study of the lateral automatic landing control laws for the Augmentor 
wing airplane is described in Reference 3. Several control law alternatives were 
evaluated in that work. Some of these results are applicable to the Twin Otter 
because both are STOL aircraft flving a 6* to 7.5* final approach at about 70 knots. 
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Differences between the powered 11ft and light wing loading aircraft In pitch 
are very significant but In roll and yaw both vehicles use the same basic control 
techniques. 

The runway alignment results that were obtained for the Augmentor Wing and 
applied to the Twin Otter are listed here and the tradeoff studies leading to 
these results are described In detail In Reference 3. 

Runway Alignment Configuration 

A forward slip technique was traded against the flat decrab, resulting In a 
clear superiority of the forward slip since the alignment Initiation timing Is 
not critical and touchdown performance, Is less dependent on vehicle and wind 
variations. 

Align Entry 

Two schemes were evaluated. One relied on a conmand rate limiter to slow 
the forward slip entry rate while permitting a reasonably high heading loop gain. 
Although landing performance with this configuration was satisfactory It had some 
drawbacks which Included overshoots of both the reference heading and the required 
steady state bank angle. Also, full crab angle was reached at a relatively high 
altitude, requiring the full wing down compensation for up to 10 seconds wnich 
may negatively impact ride quality. The altitude scheduled reference heading error 
trajectory described In Section 4 was selected as the preferred configuration as 
it eliminated the drawbacks of the rate limited align entry. 

Rudder Bias Generation 


A sideslip angle is generated while aligning the airplane with the runway in 
the presence of a crosswind. This sideslip must be held with the rudder. Three 
methods for generating the steady state rudder deflection were evaluated. An 
open loop rudder predict term, based on heading error at align entry, was evaluated 
in conjunction with the flat decrab configuration and was found to be subject to 
limitations inherent in most open loop compensation. In the second technique the 
steady state rudder command was obtained from a crosswind estimate . Good performance 
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was obtained with this method for the Augmentor Wing airplane, but the computation 
Is fairly complex and It requires a prior knowledge of several of the airplane's 
stability derivatives. The third method that was Implemented In the recommended 
control law Is a closed loop rudder control using heading reference trajectorv 
error and Its Integral to drive the rudder (as described In 4.2). The aircraft Is 
able to track the model trajectory with only small errors, thus allowing the use of 
a high gain Integral term to provide the required steadv rudder during the 
alignment. 

Wing Down Compensation 

A roll command Is needed at align entry to minimize the roll transient due to 
the align rudder kick. Three techniques were considered for the Augmentor Wing 
airplane. An open loop wing down predictor, a crosstrack acceleration referenced 
signal and one based on pseudo crosstrack acceleratirn computed from bank angle and 
rudder deflection. The last two methods produced good results and the 
accelerometer referenced signal was selected for the recommended signal because the 
Augmentor Wing airplane (and the Twin Otter) have three axis accelerometers and the 
accelerations are resolved Into runway axes as oart of the navigation 
computations. An additional compensation signal, based on rudder and yaw rate to 
roll wheel crossfeeds, was needed for the Twin Otter, as described In a following 
section. 

The following sections described tradeoff studies and alternative methods 
that were used to solve problems that were peculiar to the Twin Otter. 

6.3.2 ROLL WHEEL CROSSFEEDS 


The reconnended Twin Otter lateral /directional control law Includes 
crossfeeds of rudder coninand and yaw rate to the wheel, as shown In Figure 4-7. 
These crossfeeds were needed to minimize a pronounced tendency of the airplane to 
deviate from the localizer, away from the wind, during runway alignment. A 
sequence of alignment time histories without the crossfeeds and with 15 knots 
steady crosswind Is shown In Figure 6-4. The first time history Is with all the 
nominal Twim Otter stability derivatives. During the first two or three seconds 
Into the alignment the airplane banks In the wrong direction and a lateral devia- 
tion developes and reaches -5.18 m (-17 ft, away from the wind) at touchdown. In 
the following time histories one airplane stability derivative is set to zero for 
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Figure C>-4. The Effect of Stability Derivatives on Literal Deviation 








12 


each run in order to identify the airframe parameter that contributes the most 
to the problem. It is Quite obvious that the dihedral effect, L. , is the 

D 

most significant single contributor. 

Crossfeeding yaw rate to the wheel proved to be effective in controlling 
the adverse roll in alignment but it did not reduce the problem to a satis- 
factory low level and therefore additional crossfeeds were considered. 


Sideslip to Wheel Crossfeed 


All wheel and rudder control laws for this configuration are as shown in 
Figure 4-7 with the exception of the substitution of sideslip angle for rudder 
command, as shown in Figure 6-5. 



W ■ *•» d^c • '‘sCF ■ 0-6 -if- 

FIGURE 6-5. SIDESLIP TO WHEEL CROSSFEED BLOCK DIAGRAM 

An alignment time history with these crossfeeds is shown in Figure 6-6 and it is 
obvious that the results are very good as the lateral deviation is kept below 
0.76 m (2.5 ft). 


Sideslip angle, however, is not directly available on the Twin Otter. It 
can be computed through the use of the following equation: 


I U, 


g« 


- r dt 
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But the signal obtained this way Is susceptible to drift problems and It was 
therefore not adopted in the recommended control law. 

Heading Error to Wheel Crossfeed 

Heading change during alignment is a reasonably good aoproxlmati on for 
slldsllo In the case of a steady crosswind. An attempt was made to use 
- Atl) (see Figure 4-7) instead of 8 In the block diagram of Floure 6-F. 
Alionment time histories with 8 and used as crossfeed to the 

wheel, are given In Figure 6-7. The headlno error crossfeed results In centerline 
tracking accuracy with a steady crosswind that is as good as with the sideslio 
crossfeed. Significant excursions occur, howeyer, with a shearing crosswind and 
therefore this crossfeed has been discarded. 

Rudder Command to Wheel Crossfeed 


A crossfeed of rudder command, along with vaw rate, to the wheel produces 
oood tracking accuracy as seen In Figure 6-6 and these crossfeeds were included in 
the final configuration as shown in Figure 4-7. 

One problem associated with the use of rudder command to roll crossfeed is a 
non zero rudder trim value which could introduce a disturbance into the roll axis 
during glide slope capture. An attempt was made to incorporate a washout on 
to but even unreasonably long washout time constants (100 
seconds) resulted In a pronounced deoradation of runway alionment performance, as 
can be seen in Figure 6-8, precluding the use of a washout. The impact of rudder 
trim on landing performance is therefore evaluated. Assume that most of the rudder 
trim change on final approach is a result of thrust assvmmetrv that occurs while 
thrust is reduced in the glide slope capture maneuver. The thrust variation 
between straight flight and a 7.5® descent is: 

AT a W sin Y 

W » 5000 kg (11,000 lb) Y » -7.5® 

result in AT « -652.6 kg (-1436 Ibl 

Further assuming 10 percent thrust asswnnetrv and with the moment arm of 2.56 m 
(8.d ft) the total vawing moment is: 
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Figure 6-7. Aiignment Time Histories with Various Crossfeeds 
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N • AT * Engine • 65.3 * 2.56 - 167 kg - m (1206 ft-lb) 

The rudder deflection required to comoensate for this woroent is comouted from 


AN ♦ qsb . 0 

The result is 1® of rudder. This effect was simulated by inserting a ramp of 
rudder command in the crossfeed oath going from zero to 1.5® (for conservatism) in 
10 seconds. The lateral touchdown deviation induced bv the rudder trim change was 
determined to be less than 1.07 m (3.5 ft) for a glide slope capture maneuver at 
274 m (900 ft) height, as seen in Figure 6-9. 

6.3.3 RUDDER AND ALIGN LIMITS 


The rudder authority limit and the yaw rate command limit (Rljm in Figure 
4-7) have a significant impact on runway alignment performance. Initially, 
software rudder limits in the Twin Otter automatic landing system were set at *5® 
(out of a mechanical travel of 21® to the right and 17® to the left). Statistical 
landing results were obtained for several levels of rudder and yaw rate command 
limits and the results are summarized in Table 6-IV. Data were taken with the 
recoflWl?S3<Kl alignment control law but with (the alignment model value at low 
altitude; see Figure 4-7) of 2® and therefore a 1® mean touchdown heading deviation 
is obtained even without limits. The ±2o heading dispersion increases bv a 
factor of 4 when the vaw rate command is limited to i3.125® Instead of being 
unlimited. No degradation is apparent when the unlimited rudder authority limit is 
narrowed to ±12® and the align limit to ±2.5®/s. Table 6-V summarizes the results 
with a fixed rudder limit of ±7® and an alignment limit of il.25®/sec as a function 
of wind and turbulence level. Even with 5 knots of crosswind the touchdown heading 
dispersion is much worse than with the unlimited svstem with 15 knots of cross- 
wind. The results clearly indicate that t5® of rudder authority is inadequate. A 
12® to 15® authority is needed with a high align rate limit. The recommended 
control law uses 12® rudder authority and a 4®/sec align limit. Since the rudder 
authority is a function of dynamic pressure, data were taken for both 12® authority 
and lO® which is considered to be the worst case. The results are sunmarized in 
Figure 6-10 and the 10® author^-.' is marginal with 15 knots crosswinds. 
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TABLE 6-IV ALIGN PERFORMANCE WITH VARYING RUDDER AUTHORITY 


Align Limit Rj^jm (deg/sec) 

1 

1.25 

2.5 

3.125 

3.125 

Rudder Limit (deg) 

5 

10 

12 

15.0 

09 09 

tjQ ("leg) 


9.4 ±2.6 

4+2.5 

-3.2 

-1.1 

1.9*j-j 2±0, 


NOTE: 1. Data obtained with 15 knots shearing crosswind and corresponding 

turbulence. 

2. Results are given as u ± 2a. 


3. 


TABLE 6-V ALIGN PERFORMANCE WITH 7° RUDDER AUTHORITY 


Cross Wind (knots) 
' 1^0 (deg) 


10 


15 


1.5 5 ±2.8 9.4 ± 2.6 


NOTE: 1. Rudder align and position Hit Its of 1.25®/sec and 7® respectively 

were used. 


2. The crosswind level also defines the appropriate turbulence and 
shear. 


3. ^^-2®. 


6-24 
















The steady state rudder deflection required to handle a 15 knot crosswind Is 
readily comouted front: 

NbB + N^rSr « 0 


The result Is 11.7® which suggests that the recontmended 12® rudder limit Is the 
minimum allowable value and any lower value will be detrimental to oerformance. 

6.3.4 ADDITIONAL VARIATIONS 


Minimum Heading Reference 

The recommended alignment control law utilizes a reference trajectory for 
runway alignment as shown In Figure 4-7. This trajectory starts out with 
the ore-allgn heading error, at a gear height of 46.7? m 050 ftl and It decreases 
and levels off at the final value of at 15.24 m ^50 f tl . A non zero value 
of results In an alignment command that Is degrees shv- the runway 
heading such that the alignment Is not comolete at 15.24 m i50 ftl hut the 
Integrator reduces the heading error to zero at touchdown. This Is done to 
eliminate possible oyershoots of the runwav heading. Alignment time histories for 
the Twin Otter with 15 knot shearing crosswind are shown In Figure 6-11 with 
of 2® and 0. No oyershoot tendency Is apparent and the alrolane touches down about 
1.5® short of the runway heading with the non zero Therefore, is 
zero for the reccwnended system.* 

MLS Receiving Antenna Location 

Most of the work In this study was done assuming that the nose located 
antenna is C.G. corrected such that the input to the lateral navigation filter is 
lateral deviation at the airplane's center of gravity. Landing time histories with 
a nose mounted uncorrected antenna are given In Figure 6-12. The first three of 
the six traces are with 15 knots shearing crosswind and the second group of three 
traces Is with 15 knots steady crosswind. The first and fourth traces are with 10 
knots shearing tailwind, the second and fifth are with zero headwind, and the third 
ana sixth are with 25 knots shearing headwind. With the uncorrected antenna, the 
airplane's mse tracks the localizer and since during alignment it tends to rotate 
about its center of gravity It is displaced from the centerline at touchdown bv uo 
to 1.52 m (5 ft) away from the wind. Lateral disolacement with steady crosswind is 
larger than with shearing crosswind. 
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Increased Qalns 


The Impact of Increased gains on lateral landing deviations had been studied 
prior to the determination to use yaw to roll crossfeeds to Improve control. 

Increased from 4 to 24 deg/deq 

Kp Increased from 2.14 to 3.0 deg/deg/sec 

Increased from 0.1^7 fO.06) to 0.P56 fO.30) deq/m (deq/ft) 

Time histories with the Increased kr are qlven In Figure 6-13. 

These gains Improve lateral control but a touchdown deviation of 2.13 m (7 ftl 
still occurs with 15 knot shearing crosswind. Time histories with the Increased 

Kr and Ky are given In Figure 6-14. Control Is Improved but It Is not 
so good as with the cro^sfeeds (see Figure 6-1). Also, these high gains result In 
excessive wheel activity and they are therefore not utilized In the final svstem 
configuratl on. 

The label "Hold at T.D." In Figure 6-13 means that the simulation computer Is 
switched to the Hold Mode at touchdown which is the way the simulation is normally 
ooerated. "No Hold at T.D." means that the simulation was allowed to run oast the 
touchdown altitude as If the runway was not there, In order to evaluate system 
dynamics by allowing more time In which transients can develop. 
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7.0 CON g US IONS 


The conclusions of this light wing loading STg aircraft autanatic landing 

study results are given here. These conclusions are based on results obtained from 

simulation. * 

1. This light wing loading airplane ooerates on the front side of the required 
power curve and therefore a conventional control strategy can be used. 

Flight path Is effectively controlled with the elevator and airspeed with the 
throttles. This Is In contrast with the situation of the powered lift STOL 
airplane operating on the back side of the power curve. There, the throttles 
are used for flight path control and the elevator for airspeed as described 
in Reference 1. 

2 . The control laws that were developed are capable of producing landing 
performance that Is consistent with Category III operation Into STg runways 
as defined in Reference 

3. The use of the fast responding DLC spoilers enhances flight path control 
bandwidth and imoroves landing performance. The Improvement Is mostly in 
sink rate for the constant flare height configuration and mostly In range for 
the variable flare height control law. 

4. All the control laws that were evaluated resulted in a flat touchdown 
attitude with this light wing loading airplane. Attempts to arbitrarily 
command a higher attitude resulted in excessively long landings. The lift 
reduction capability of the Twin Otter spoilers was aoproximatelv 0.13 g. A 
higher lift reduction capability, if available, would have probably allowed 
rotation to higher attitudes without floating. Lower aooroach airspeed for 
spoilers might have also been beneficial. 

5. The constant flare height closed loop longitudinal control laws produce equal 
or better landing accuracy and are much simpler in structure than the 
variable flare height control laws, using predictive oitch and elevator terms. 

K. Landing performance results that were obtained for the light wing loading 

Twin Otter are very similar to these obtained in a previous studv (Reference 
11 for the Dowered lift Augmentor Wing airplane despite almost a factor of 
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two difference In wing loading between the two aircraft. The two aircraft 
approach to land at approximately the same speed. 

7. Turbulence and winds are the major contributors to landing dispersions. 
Beam disturbances have no significant Impact on the longitudinal axis 
but they do affect lateral performance. 

8. Yaw rate and rudder to wheel crossfeeds were used to overcome the adverse 
effect of a pronounced dihedral effect on lateral landing performance. 


7-2 


APPENDIX A 


SIMUUTION DEFINITION 


TABLE OF CONTENTS 
APPENDIX A 


Airframe Dynamics 

Page 

A-1 

A-1.1 

Equations of Motion 

A-2 

A-1.2 

Stability Derivatives 

A-2 

A-1.3 

Airframe Response Characteristics 

A-2 

Control 

System Dynamics 

A-2 

Geometry and Sensors 

A-12 

A-3.1 

Sensor Geometry 

A- 12 

A-3.2 

Physical Data 

A-15 

A-3.3 

Sensor Models 

A-15 

Beam Disturbance Models 

A-17 

A-4.1 

Beam Noise 

A-17 

A-4.2 

Azimuth Signal D1*scret1zat1on 

A-19 

Atmospheric Disturbance Models 

A-19 

A-5.1 

Standard Wind Model 

A-21 

A-5.2 

Logarithmic Wind Shear 

A-23 

A-5.3 

Altitude Profiles for Lateral Landings 

A-23 



I 



LIST OF FIGURES 



APPENDIX A 


FIGURE 

TITLE 

PAGE 

A-IA 

Free Airframe Responses to Control Steps - 
Longitudinal 

A-7 

A- IB 

Free Airframe Responses to Control Steps - 
Longitudinal 

A-8 

A-2 

Free Airframe Responses to Hind Steps - 
Longitudinal 

A-9 

A-3A 

Free Airframe Responses - Lateral /Directional 

A-10 

A-3B 

Free Airframe Responses - Lateral/Directional 

A-11 

A-4 

Longitudinal Actuator Models 

A- 13 

A- 5 

Lateral -Directional Actuator Models 

A- 14 

A-6 

Airplane Geometry 

A-15 

A- 7 

MODILS and MLS Error Models 

A-18 

A-8 

MODILS Signal Discretization 

A-20 

A-9 

Standard Wind Model 

A-22 

A- 10 

Lii<ear and Logarithmic Shear Profiles 

A-24 

A-11 

Altitude Profile Generation 

A-25 

A- 12 

Typical Altitude Profiles 

A-26 


1 

i 


1 


A - i1 


LIST OF TABLES 
APPENDIX A 


TABLE TITLE PAGE 

A- I Longitudinal Equations of Motion A-3 . 

A-II Longitudinal Stability Derivatives A-4 

A- I II Lateral /Directional Equations of Motion A-5 

A-IV Lateral Directional Stabil ity Derivatives A-6 

A-V Physical Data for Landing Configuration A-16 

A-VI Sensor Characteristics A-17 


A -iii 


APPENDIX A 


SIMUUTION DEFINITION 


The non-plloted simulation of the modified DeHavilland DHC-6-100 Twin 
Otter research STOL aircraft and its automatic landing system are defined 
in this appendix. This simulation definifon includes: 

1. Airframe Dynamics 

2. Control Dynamics 

3. Geometry and Sensors 

4. Beam Noise Model 

5. Wind Models 

A-1 Airframe Dynamics 


The three nominal 

flight conditions for 

this study are 

defined below. 


LONGITUDINAL 

LATERAL/DIRECTIONAL 


NO DLC 

DLC 


Weight, kg (lb) 

4990 (11000) 

4990 (11000) 

4990 (11000) 

Speed, kt 

71 

75 

70.5 

Angle of Attack, deg 

-3 

-3 

-3.8 

Glide Slope, deg 

-6 

-6 

-7.5 

Flaps, deg 

35.4 

35.4 

35.4 

Spoilers, deg 

0 

20 

0 


It should be noted that the glide slope angle for the longitudinal studies 
is 6“ whereas 7.5" was used for the lateral/directional work. This is due 
to the fact that the lateral studies were conducted earlier when 7.5" was 
planned to be used as the nominal approach angle. Later on, the 6" glide 
slope was selected as a more appropriate value. 

The approach speed with spoilers is higher in order to provide adequate 
stall margin with the wing spoilers extended. 
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A“l.l Equations of ftotlon 


The nonnal set of uncoupled, linearized, small perturbation aero- 
dynamic equations of motion were used as documented in Tables A- I and A-III. 
The aircraft equations and stability derivatives included in the separate 
longitudinal and lateral landing simulations are described in detail along 
with representative time responses. Note that the important aerodynamic 
nonlinearities, such as ground effect, were also included in the simulation. 

A- 1.2 Stability Derivatives 


Twin Otter longitudinal dimensional stability derivatives and other 
pertinent data are included in Table A-II and the lateral directional data 
are in Table A-IV. The aerodynamic data are based on NASA's simulation 
documentation. Spoiler derivatives were obtained from flight tests by NASA. 

A-1.3 Airframe Response Characteristics 

Characteristic roots for the nominal approach cases are given below. 


Longitudinal (75 kt) 

Sp“ 0.605 

“sp‘ 

2.24 rps 


Sh* 

“ph* 

0.313 rps 

Lateral (70.5 kt) 


“or* 

1.467 rps 


■ -11.27 sec 

^R ’ 

0.246 sec 


The longitudinal free airframe responses to step elevator, throttle and DLC 
spoilers are given in Figure A-1 and responses to u and a gust inputs are 
given in Figure A-2. Lateral /directional free airframe responses to wheel, 
rudder and 6 gust steps, as well as roll rate and yaw rate initial condition, 
are given in Figure A-3. 

A-2 Control System Dynamics 


For the Twin Otter vehicle, control surface aerodynamic and inertia loaus 
were sufficiently small that acceleration limits and detailed actuator models 
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TABLE A-I. LONGITUDINAL EQUATIONS OF MOTION 


a » 

w » 

q = 

*NACC * 

9 ^ 

*XACC 


*A '*co*yo HF + Xg^ 6 ^ +X5SP 6sp 
“a Z;, ^^ + Zw '^^A ■*'(U, + Z^)q+AZQg HF + Z5, 6, +Zs^ 6* + Z«sp «sp 

q * AMgg HF + Mj, 6, +Mg^ 6^ 6sp 

(U,q-w)cosofQ-ilsinofQ q 

u + g(# 


•’CG * U, (0^ -OTo)-W, COSTq 
^G * hcQ “ ^ ^G ®T . 

*^RA “ *^G 

^8 “ *^G ’ Vo ^G “>0 > “ 2 rG =°S ^0 


“a “ ^1 “ ^O “wWD 

'*'a = ^+«^WIND 

HF = c'^'g ^*^ge 

NOTES; are incremental aerodynamic values about trim. 

(f is incremental pitch attitude about trim. 

the subscripts “o ” and t indicate trim and total values respectively. The subscript i indicates 
inertial quantity 
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TABLE A-II. LONGITUDINAL STABILITY DERIVATIVES 


W kg, (lb) • 4990 (11,000), • 

ly kg-m^ (slug-ft^) - 33082 (24400) 


- 35. 4* 


WITHOUT SPOILERS 


WITH SPOILERS 


kt 

deg 

1/sec 

1/sec 

2 2 

m/sec /rad (fps /rad) 

2 2 

m/sec /rad (fps /rad) 


71 

0 

-0.0883 

0.1346 

4.356 (14.29) 


75 

20 

-0.0911 

0.1423 

4,866 (15.96) 
-0,3573 (-1.172) 


1/sec 

1/sec 


2 2 

m/sec /rps (fps /rps) 

2 2 

m/sec /rad (fps /rad) 

2 2 

m/sec /rad (fps /rad) 

m/sec^/rad ( fps^/rad) 


-0.4836 

-1.009 

-0.0086 

-0.953 (-3.127) 
-2.896 (-9.50) 
-5.360 (-17.59) 


-0.4578 

-1.066 

-0.0086 

•1.007 (-3.304) 
-3.234 (-10.61) 
-5.989 (-19.65) 
3.566 (11.70) 


M* 

u 

M 

2 2 
rps /m/sec (rps /fps) 

0.00478 (0.00146) 

0.00533 (0.00163) 

2 2 * 
rps /m/sec (rps /fps) 

-0.09161 (-0.02792) 

-0.09677 (-0.02949) 

w 

M- 

w 

rps^/m/sec^ (rps^/fps^) 

-0.1006 (-0.00306) 

-0.1006 (-0.00306) 

M 

1/sec 

-1.203 

-1.271 

q 

M. 

1/sec^ 

-3.384 

-3.778 

6e 

M. . 

1/sec 

-0.6096 

-0.6808 

6th ^ 


1/sec 

- 

0.213 

6Sp 


2 2 
m/sec (fps ) 

0.179 (0,589) 

0.200 (0.657) 


m/sec^ (fps^) 

•y 

-1.02 (-3.36) 

-1.14 (-3.75) 

^GE 

rps 

-0.551 

-0.615 

"ge 

m (ft) 

3.66 (12) 

3.66 (12) 


NOTE: ALL DERIVATIVES ARE IN STABILITY AXES 
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TABLE A-III. LATERAL/DIRECTIONAL EQUATIONS OF MOTION 


Lateral /D1 recti onal 








"^A • 'a 


P * '■d * ^A ■*■ + L^ ♦ r, + L 

^ X 


r ' ' A ■ ‘■0 • 6 *’8*®A 


'■‘-6R * ^R ^ '-SA * ^A 


I 




+«S-®*"6-®A 


^6R * ^6A • 


PI i COS ajj - Sin aQ-j | p 


r B sin cos a 


oJ 


Body axes rates 


= Pg + tan 0Q . rg 


Euler rates 


*y * ^0 '"b) ■ 9*^ 

Yr * (dy + 9(J|) COS ♦ cos Av 
®A “ ®WIND 


Lateral acceleration 


Runway crosstrack acceleration 


Ai|; s i|; - 4/ 


RUNWAY 


The A and o subscripts denote aerodynamic and trim quantities, respectively. 
The B subscript denotes body axes quantities. 

The subscript I indicates inertial quantity. 


TABLE A-IV. LATERAL DIRECTIONAL STABILITY DERIVATIVES 


% 

kt 

70.5 

“o 

deg 

-3.8 


kg-m^ (slug-ft^) 

22776 (16800) 

iz 

kg-m^ (slug-ft^) 

50908 (37550) 

^xz 

kg-m^ (slug-ft^) 

1898 (1400) 


- 

-0.00473 

^r/^o 

- 

0.0236 


1/sec 

-0.1756 


j 1/sec 

0.068 

>6a/^c 

j 1/sec 

0.00059 

Lp 

1/sec 

-4.089 

Lr 

1/sec 

2.9896 

l: 

1/sec 

0 

B 

■■sR 

1/sec^ 

1/sec^ 

-2.052 

1.321 

•-6A 

1/sec^ 

-5.473 

% 

1/sec 

0.0067 

Nr 

1/sec 

-0.6880 

Nb 

1/sec 

0 

Nb 

1/sec^ 

1.75 

'‘aR 

1/sec^ 

-1.818 


1/sec'^ 

-0.4872 


NOTE: 1. All derivatives and inertias are given in stability axes. 

2. All angles are in radians unless otherwise noted. 
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Figure A-3B. Free Airframe Responses • Lateral/Directionsl 



were not re<)u1red. Thus only Important non! inear1t1es as rate and position 
limits were Included In this study. The simplified actuator models are 
described In this section. 

The throttle, engine, elevator and spoiler actuator models are given In 
Figure A-4. The one second lag In the throttle servo, as shown In the Figure, 
was Included with the variable flare height and the low gain constant flare 
height control law configurations. The high gain constant flare height con- 
figuration assumed a high bandwidth throttle servo btc the same engine dynamics 
as shown In the Figure. The turboprop engine dynamics are different for 
throttle advance or retard. The throttle and elevator position limits of 
Figure A-4 are with respect to a zero wind nominal trim position. Spoiler 
limits are In terms of actual sp< 'ler deflection. 

The wheel-aileron and rudder models are shown In Figure A-5. Note that 
Initially mechanized rudder limits had been ±5° and were found to be Inadequate 
for the 15 knot crosswind landing requirement. Accordingly, the rudder limits 
were modified to The limits are dependent on dynamic pressure and a ±10° 

worst case limit was assumed during part of the simulatfSflf J 

A-3 Geometry and Sensors 

A-3.1 Sensor Geometry 


The relative geometry of the gear, c.g., and MLS antenna are illustrated 
in Figure A-6. The gear and MLS receiver location are expressed in terms of 
eg height above the runway by the expressions: 

hg • h^g - Zg cos 9 Xg sin 6 ■ 5.97m il9.58 ft) 

^REC " ^cg ■ %C ® * ^REC ® ^REC * 


A-12 
















FIGURE A-6. AIRPLANE GEOMETRY 


A-3.2 Physical Data 

Sensor location, landing gear geometry and the weight, moments of inertia 
and dimensions of the Twin Otter airplane in the landing configuration are given 
in Table A-V. This geometry defines the following absolute touchdown constraints 
for this vehicle: 


e-|-g max 

120 

min 

-1.5 

wing scrape 

19“ 


A- 3. 3 Sensor Models 


Only those sensors whose dynamics or errors impact landing performance 
are discussed in this section. Also, only the errors which affect landing 
performance are included. The properties of the sensors which impact landing 
performance are summarized in Table A-VI. 


Although MLS yields discrete information at 5 scans per second for 
elevation 1 and azimuth guidance, and 40 per second for DME, continuous pos- 
ition inputs were used during these studies. During the previous glide slope 
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TABLE A-V. PHYSICAL DATA FOR LANDING CONFIGURATION 



Sensor Location 

11 

WL 

BL 


CG (26X MAC) 

529.64 (208.52) 

254 (100) 

0 


MLS Antenna 

-67.31 (-26.5) 

177.8 (70) 

0 


Radar Altimeter 

589.28 (232) 

127 (50) 

0 


Accelerometers 

458.88 (180.66) 

254 (100) 

0 

GEAR GEOMETRY 

Nose 

Right 


Left 

^G 

3.937 (12.9) 

-0.596 (-1.956) 

-0.596 (-1.956) 

^G 

0.0 

1.828 (6.0) 


-1.828 (-6.0) 

Zg 

1.99 (6.53) 

1.99 (6.53) 


1.99 (6.53) 

^Compression 

.253 (0.83) 

.311 (1.02) 


.311 (1.02) 


AIRCRAFT DIMENSIONS 



W 

48947 N 

(11,000 

1b) 


22776 Kg-m^ 

(16,800 

s1ug-ft^) 

^Y 

33079 Kg-m^ 

(24,400 

s1ug-ft^) 

iz 

50908 Kg-m^ 

(37,550 

s1ug-ft^) 

^XZ 

1898 Kg-m^ 

(1,400 

s1ug-ft^) 

^wing 

39.02m2 

(420 ft2) 

b 

19.812 m 

(65 ft) 

r 

1.98 m 

(6.5 ft) 


NOTES: 1) Sensor location 1s expressed In centimeters (Inches) 

2) Gear geometry Is expressed In meters (feet) 

3) WLCG Is estimated. 


and localizer track MLS studies reported In Reference A-1, It was determined that 
these update rates provided control activity and landing performance Identical to 
continuous guidance signal, especially If beam filtering Is used. Although 
continuous position signals were used to limit simulation complexity, the actual 
MLS error models defined below were Included to maintain fidelity In the results. 
In addition, the effect of guidance signal resolution was Investigated using the 
model described In the following section. Only azimuth effects were considered 
since the effects of the elevation siqnal resolution are minimal. 

TABLE A-VI. SENSOR CHARACTERISTICS 


SEflSOR 

DYNAMICS 

ERRORS 

Radar Altimeter 

1 

71TTT 

Bias - ±1.03m f±3.4 ft) 

GS/IOC Receiver 

.15 + 1 

Figure A-7 

Vertical Gyro 


Verticallty - ±.6® 
False erection - ±1.0® 

Accelerometers 


Cross axis sensitivity .01 

Course Datum 


Equivalent Bias - ±4.0® 


NOTE: All errors are given as 4.5a values. 
A-4 Beam Disturbance Models 


1-4.1 Beam Noise 


In this stuty, both MOOILS and HLS error models were used. The azimuth, 
elevation, and DME error amplitudes, spectral characteristics, and resolution are 
qlven in Figure A-7. It should be noted that the effect of DME inaccuracies is 
small with respect to the anaular errors durino the final approach. The noise and 
bias values for the MLS are based on the preliminarv specification for the system 
at the Navv Crows Landing Auxiliary Landing Field. 


ELEVATION 


'EL 
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1 

ERROR (RAD) 

V 

NOISE 

telS+1 


A 




RANGE 


AZIMUTH 


WHITE 

NOISE 


'A2 



>AY 
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‘^NOISE 

*^BIAS 

RESOLUTION 




MODILS 



Flevation 

O 

• 

o 

sec 

0.07* 

0.05* 

0.01* 

Azijnuth 

4.0 

sec 

0.03* 

0.17* 

O.l* 

DME 

- 


12.19m (40 ft) 
• 

6.10m (20 ft) 

18.29m (60 




MLS 



Elevation 

0.5 

sec 

0.0625® 

0.0375* 

O.Ol® 

Azimuth 

0.5 

sec 

0.105* 

0.0935* 

O.Ol* 

DME 

0.5 

sec 

12.19m(40 ft) 

15.24m (50 ft) 

1.83m (6 ft) 


FIGURE A-7 MOOILS AMD MLS ERROR MODELS 




A-4.2 Azimuth Signal Discretization 


Two significant sources of discretization exist In the lateral guid- 
ance signal : 

1} Azimuth resolution of .1 degrees for MODILS, .01 degrees for MLS. 

11) Digital autopilot scaling of 1.22 m (4 ft) per bit. 

The MODILS lateral deviation resolution Is a function of slant range. 
Assuming the typical STOL-port geometry shown In Section 3, with a 7.5 degree 
descent path, the .1 degree discretization yields about 1.07 m (3.5 ft) at 
touchdown and 6.1 m (20 ft) at 304.8 m (1000 ft) altitude. Thus the azimuth 
resolution Is the major contributor to lateral discretization with MODILS, 
and Its Impact on track and align performance was determined. The Implementa- 
tion diagram and discretization effects with sinusoidal Input are given as 
Figure A-8. 

With MLS, the azimuth resolution Is Improved by an order of magnitude. To 
take advantage of this Improved guidance signal quality, the digital signal 
resolution was Increased for approach, A resolution of better than .305 m 
(1 ft) would provide an essentially continuous guidance signal. Thus a con- 
tinuous azimuth signal adequately models the MLS system characteristics. 

A DME resolution of 18.532 m (.01 N. miles) yields only .305 m (1 ft) 
discretization error at a 1 degree azimuth error. Since aircraft azimuth 
deviation Is always less than 1 degree during final track and alignment, the 
DME resolution has insignificant impact on lateral signal discretization. 

A-5 Atmospheric Disturbance Models 

In these landing studies, both a standard atmospheric disturbance model 
and specific deterministic wind Inputs were considered. The final performance 
determination was based on the standard disturbance model. 
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A-5.1 Standard Wind Model 


The wind model used for obtaining all the statistical data Is closely 
patterned after the standard FAA wind model specified In AC-20-57A (Reference 7) 
and described more fully In References A and 5. The total wind level also 
determines the turbulence amplitudes, while the shear corresponds to the headwind 
and crosswind components. A summary of the standard wind model Is shown as Figure 
A-9. The wind shear of this model consists of a steady wind at altitude down to 

61 m (200 ft.) at which point the wind decreases linearly with decreasing altitude 

all the way to the ground. The nominal wind magnitude is given as the value at the 
7.62 m (25 ft) reference altitude. Thus, a 25 knot shearing headwind has a 42.5 

knot magnitude at altitude, decreasing to 22.5 knots at zero height. Similarly, a 

10 knot shearing tailwind decreases In magnitude from 17 to 9 knots and a 15 knot 
crosswind shears from 25.5 to 13.5 knots. The magnitude of the horizontal 
turbulence, according to this model, is proportional to the wind level, as shown In 
Figure A-9. For conservatism, most statistical data In this study were taken 
assimlng a 70 percent orobabllltv of encountering a 25 knot shearing headwind with 

the associated of 3.75 knots and a 30 percent orobabllltv of encountering a 

10 knot shearing tailwind with a turbulence level of 2.70 knots RMS. The 

horizontal turbulence levels given above were computed from the formulae m 

u 

0.15 WINOV given in Figure A-9. WINOV 1s the total wind velocity and In the case 
of the 25 knot headwind, WINOV Is 25 knots. The 10 knot tailwind Is assimed to be 
associated with a 15 knot crosswind, resulting in a total wind of 18 knots. These 
maximum conditions are referred tt) as limiting winds In the body of the report. 

The localizer track and runway alignment data were taken with 15 knots shearing 
crosswind and of 2.25 knots. Reference 7 assigns a 1 oercent Drobability 
with headwinds In excess of 25 knots and 4 percent with crosswinds exceeding 15 
knots. The much higher probabilities used in this study produce conservative 
results. A constant level of 1.5 knots RMS vertical turbulence (invariable with 
altitude or wind conditions) was used. For this study, winds were assuned to be in 
earth local level axes and transformed into aircraft axes. No pitch rate gusts 
were used, since their effects are negligibly small comoared to horizontal and 
vertical turbulence. Uncorrel ated white noise generators were used for 
longitudinal, lateral, vertical, and roll rate gusts. 


A-?l 


ia- 


u> 



^ * n ®Q [•“ 




1 

1 


WHITE N. 

\ 

Ji J 

1 


NOISE ^ 

Tx, 


'SPAN+1 








U, 


wind 




cot^i - lin^i 






WHITE N. 

1 

^ - 

«w 

k u 

1 


NOISE ^ 



’^AREA®'*^ ^ 



1 ■»^aoo4 


<h*CQ-W 


P- 





!v_ 

WHITE -v 

Uo 

L_J 

1 X 
1 ^ 


NOISE ■VS' 


n 

n 






WHITE 

NOISE 


> 



.W. 


wind 


V . 
fc wind 


gutt 


gun 



u 

V 

W 

P 

•-(SCALE LENGTH) 

182.80m (600 ft) 

182.88m (600 ft) 

9.14m (30 ft) 


’’SPAN. AREA 

27.66Km 


12J9K 

m 




•Jo 


0 

aiswiNov 

ai8V 

a773^ '••Skt) 

00063— 5^^ 

•^ »c 

Mfian Wind Limit 

>1X88 >26kt 

t7.73m/iae(±16kt) 


*S.1Sm/»c -10 M 

1% 



^rob of ExoQodmoi 

4.5% 




I 


'm 


b 

U« 


CG (200 ft mwimuni) 

Wind ftofsrvnot Altitiidi * 7.62m (2B ft) 

Mnric Constant ' - 0.304tm (1.0 ft) for mathc (Cnglhh) untti 

W)n| Span “ 1fil8lm (06 ft) 

Approacf) 4>M<t 

OoM«n Wind Spaad; V “ Cro t aind S pai d ; WINDV ■ Total Wind Spaad (25 kt Limit) 


Figure A-9. Standard Wind Model 













A-5.2 Logarithmic Wind Shear 




The performance of the various glide slope track and flare control laws 
was evaluated deterministically with the standard linear wind shears described 
above and also with the logarithmic shears, as defined by the following ex- 
pression: 


“wind " “ref to -^512 logjQ (3.28084*h) + 0.3692] 
h is expressed in meters. 

7.62 m (25 ft) 


Twenty five knot standard and logarithmic shear profiles are shown in Figure A-10. 
The logarithmic shear is a more severe disturbance during the flare because of 
its steeper gradient below 15.24 m (50 ft). 

A-5.3 Altitude Profiles for Lateral Landings 

Since the time between align initiate and touchdown can have a significant 
impact on landing performance, the effect of longitudinal winds and shears on 
lateral performance was included in the simulation. 

A simplified flare model was constructed, with the altitude trajectory 
varying with inertial velocity and flare time constant in a manner very 
similar to the actual pitch approach and flare control system. The altitude 
profile generator block diagram is given as Figure A-11, with sample profiles 
for limiting headwind and tailwind shown in Figure A-12. 

This altitude trajectory is used to drive the sidewind shear and the 
align model, and to indicate touchdown. Thus the proper relationship is 
maintained between altitude and time for all downwind conditions, to allow 
realistic determination of lateral landing performance. 
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SUPPLEMENTARY GLIDE SLOPE TRACK 


AND FLARE SIMULATION RESULTS 
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APPENDIX B 

SUPPLEMENTARY GLIDE SLOPE TRACK AND FLARE 
SIMULATION RESULTS 


This Appendix contains simulation generated landing time histories and 
flare trajectories that supplement these given in Section 5 to cover all the 
major longitudinal control law configurations and deterministic disturbances 
that were evaluated. The simulation generated probability distributions of 
longitudinal landing variables, on which the Section 5 summaries are based, 
are also given here. 

The landing time histories of the variable flare height configuration and 
the low gain constant flare height configuration, both with spoilers and stan- 
dard wind shears (as defined in Appendix A) were given in Section 5 as Figures 
5-1 and 5-2. Landing time histories for all the remaining combinations of 
evaluated configurations and deterministic disturbances, including logarithmic 
wind shears (defined in Appendix A) are given here as Figures B-1 through B-10. 

Similarly, the h/h flare trajectories for the variable flare height con- 
figuration and the low gain constant flare height configuration, both without 
spoilers and with logarithmic wind shears, are given in Section 5 as Figures 
5-3 and 5-4. The flare trajectories for the remaining configurations and dis- 
turbances are given here as Figures B-11 through B-20. 

Figures B-21 through B-liO are the probability distributions of the longi- 
tudinal landing variables. These variables are touchdown sink rate, range 
(measured from the GPIP), pitch attitude and airspeed, and the deviation from 
the glide-slope at the 30.48 m (100 ft) window height is also given. In 
Figures B-21 through B-50 one curve is given per variable, obtained by com- 
bining results obtained for 25 knot shearing headwind and a horizontal turbu- 
lence level of 3.75 kt rms, with results for 10 knots shearing tailwind and 
a horizontal turbulence level of 2.70 knots. Both headwind and tailwind results 
were obtained with 1.50 knots of vertical turbulence and MLS beam noise. The 
results were combined assuming 70X probability for the headwind and 30% for 
the tailwind. Figures B-51 through B-110 have two or three curves for each 
variable, where each curve represents one specific set of deterministic and 
stochastic disturbances, as indicated on the figure or in the list of figures. 
All statistical data for the Twin Otter were obtained with the stancard wind 
models that are defined in Appendix A. 
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Figure B-89. Constant Flare Height. No Spoilers, Low Gains 
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Figure B-91. Constant Flare Height, High Gains, Spoilers 
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Figure B-93. Constant Flare Height, High Gains. Spoilers 
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APPENDIX C 


SUPPLEMENTARY LOCALIZER TRACK AND RUNWAY ALIGNMENT 
SIMULATION RESULTS 


This Appendix contains a discussion of the effects of gyro and accelero- 
meter errors on lateral landing performance. Simulation generated landing 
time histories to supplement these given in Section 6, are also included here. 
Landing time histories are given for the various control law variations that 
were evaluated during these studies. Localizer track activity traces and 
time histories of responses to various failures are given, as well as responses 
to sensor errors. Simulation generated probability distributions of the 
lateral /directional landing variabler are also given here. Probability distri- 
butions for various system and oisturbance variations are included. Perfor- 
mance results that are siminarized in Section 6, are based in part on these data. 
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C-1 DISCUSSION OF GYRO AND ACCELEROMETER ERROR EFFECTS ON LATERAL AXIS 
PERFORMANCE 


The lateral axis control law accepts as Inputs the following signals: 

1. Roll attitude 

2. Roll rate 

3. Runway axis lateral acceleration 

4. Pos and rate data from the navigation filter (Y^, Y^) 

These signals are necessarily corrupted by a finite amount of noise and 
errors of various kinds. The most commonly discussed noise Is that associated 
with the localizer deviation signal and is the reason for the use of the navi- 
gation filter. The sensitivity of the autoland system to this source of noise 
1$ assessed with the simulation by applying random gausslan noise through the 
appropriate filter and summing this with the Y deviation data. The attitude, 
attitude rate, and acceleration signals are normally characterized as havinc a 
bias (null error). Slow time varying attitude errors due to gyro drift and 
erection to a false vertical also must be considered. Furthermore, the runway 
axis lateral acceleration signal, which Is used for wing down compensation in 
align and for localizer signal augmentation in the navigation filter. Includes 
signal product terms such as a^^ Ay and therefore will contain time varying 
error signals due to bias errors In either signal. 

The following material will evaluate these errors quantitatively and a 
summary of the analysis will be offered at the end of this section. 

C-1. 1.1 ACCELEROMETER ERRORS 


Accelerometer errors effect th» aircraft through the path. The expression 
for y*p was developed in Reference 3 where it Is shown that the dominant terms are 

wnere the prime Indicates the error components and the 1 indicates the correct 
value. Each of these terms is being effected by accelerations occurring at 


C-2 


different times: 




1 


a^' - occurs mainly due to accelerometer cross axis sensitivity and there- 
fore will be affected whenever sustained forward or normal accelera- 
tion exist. Such accelerations may occur under the following conditions: 


a) During the pushover manuever (H < 1000 ft) 

Mainly the forward acceleration is affected, and the estimated 
value is a^^^ ■ 1 ft/sec^. It is also assumed that this change 
occurs as a step for a duration of 20 sec (estimated time for the 
pushover maneuver to be completed). 

b) Due to wind shear 

According to the defined wind model, wind shear will affect 

the forward acceleration starting at 200 ft through touchdown. 

The estimated value is a„ ■ 1 ft/sec^. 

’‘I 

c) During Flare (H < 50 ft) 

Ouri.ig this p*hase, the normal and longitudinal accelerations 
are present and the maximum estimated value is 

^ ^t/sec^ (4.5o ) 
a^ = 2 ft/sec^ 


a 


’‘I 


AiJ) • 


The maximum (4. So) estimated compass systen* error is 4". This 
term will affect the aircraft whenever linear accelerations 
(axj) exist according to the list above. 




The maximum esti.iwted roll angle error is 1.6° (4.5®). This term 
will affect the aircraft during the flare maneuver only because 
of the associated normal acceleration. 
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Figure C-17 shows the effect of a .5 fps^ step in error Inserted at 

an altitude of H ■ 1000'. In the first two traces, the step was rennved 

after 20 sec and in the other two the step was maintained through touchdown. 

The worst case, which is a combination of sustained acceleration error- and 

2 

shearing crosswind is 70 ft/fps and -will be assessed on a 4.Sc> basis. Figure 
C-18 shows the effect of steps inserted at altitudes of H * 200 feet (effect 

of shear) and H * 50 feet (flare) respectively. The touchdown deviation is 

2 2 

60 ft/fps and 8 ft/fps , respectively. 


Errors Due to the a 


£ 


term 


The cross axis sensitivity is estimated to be 0.01 
Therefore, the total lateral deviation will be: 


a) 

Pushover maneuver: 

1 ft/sec^ 0.01 70 » 

b) 

Wind shears: 

1 ft/sec^ * 0.01 * 60 « 

c) 

During flare: 

5 ft/sec^ * 0.01 8 » 



2 ft/sec^ * 0.01 ♦ 8 » 


RSS total 


g/g on a 4.5o basis. 


0.7 ft 

0.6 ft 

0.4 ft 
0.16 ft 
1.02 ft (4.5o) 


Errors Due to a^j a^'term 


a) Pushover maneuver: 

b) Wind shears: 


1 ft/sec^ * -g^ 3 * 70 » 4.9 ft 
1 ft/sec^ * -§^ 3 * 60 * 4.2 ft 
RSS total 6.45 ft 


Errors Due to a^j 4 "term 


This term exists only during flare; to give: 

5 ft/sec^ * * 8 • 1.12 ft 


The different contributors to the touchdown error are also RSS ed to give a 
total estimated touchdown deviation of 6.62 ft on a 4. So basis or 2.95 ft on a 
2o basis. 
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C-1.1.2. GYRO ERRORS 


Gyro errors may arise because of: 

a) Vertical ity error 

b) False erection 


Gyro data suggests max. vertical ity error of 0,6 degree (4.5a) and a 
maximum error of 1. degree due to false erection to give a total of 1.60* 
(4.5a). It is assumed that the rate of erection is in the order of 1.0 to 
2.0 degrees per minute. 


False erection may occur due to uncoordinated maneuvers which cause lateral 
accelerations, followed by erection to a false vertical. Such maneuvers may 
occur prior to final approach and during the align maneuver. The same maximum 
errors, as stated above, will be assumed for both cases. It is also assumed 
that final approach and localizer mode initiation occur at an altitude of 
H » 1000 ft. 


Figure C-19 shows the response to an existing 1.0 degree error prior to 

*C 1 

localizer mode initiation. A steady state defined by y * 


is maintained. The first trace shows the response following the mode initiation 
and the second trace shows the same, this time with the gyro erecting back at 
the rate of 1.0 degree per minute. If a worst case of tailwind is assumed, the 
elapsed time from H » 1000 ft to touchdown will be about 60 sec which will leave 
a lateral deviation error of 6 ft/degree. As a point of interest. Fig C-20 
shows the same response carried throughout the align maneuver with limiting 
crosswind and crosswind shear. The touchdown error is 5 ft/degree. 


Fig C-21 shows the same response, this time with a gyro erection rate of 
2®/min. Although the lateral deviation is peaking earlier, the same 6 ft/degree 
touchdown deviation renains after 60 sec. 


Fig C-20 also shows the response to the erection cycle occurring at the 
align altitude (H » 150 ft). It is assumed that the aircraft is centered on 
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the beam prior to align and the only disturbance Is due to the rate of erection i 

trying to respond to the lateral acceleration. The max. error at touchdown is 
about 1.5 ft. For reference, the same response with limiting crosswind and 
crosswind shear Is also shown. 
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Summary of Gyro Errors 


Errors existing prior to mode Initiation: 

6 ft/deg X 1.6* « 9.6 ft 
Errors at align: 

1.5 ft 

RSS total 9.72 ft (4.5c) 

C-1.1.3 SUW1ARY OF TOTAL ERRORS (Accelerometer and Gyro) 

Total accelerometer errors: 6.62 ft 

Total gyro errors: 9.72 ft 

RSS total: 11.76 ft 


or, a total touchdown error of 5.23 ft on a 2o basis. 


C-1.2 INCREASED Ky 


Since a higher lateral deviation proportional gain 
error analysis Is performed for * • 2Ky. Figures C-22 
most of the traces for this higher gain and the following 
tion based on the same assumptions as were defined above. 


(Ky) was considered, 
through C-24 repeat 
is the error calcula- 


C-1.2.1 ACCELEROMHER ERRORS: 


Cross Axis Sensitivity 

.) Pushover maneuver: 1 ft/sec^ X 0.01 X 48 ft/ft/sec^ * .48 ft 

(continuous step, worst case) 
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.46 ft 


b) Wind shears (H < 200 ft): 1 ft/sec^ X 0.01 X 46 ft/ft/sec^ ■ 

c) During flare: 5 ft/sec^ X 0.01 X 5 ft/ft/sec^ ■ .30 ft 

2 ft/sec^ X 0.01 X 6 ft/ft/sec^ - .12 ft 

RSS total .74 ft 


Errors Due to the aj^j Aip^term: 

1 ft/sec^ X ^ ^ ft/ft/sec^ * 3.35 ft 

1 ft/sec^ X X 46 ft/ft/sec^ - 3.21 ft 

RSS total 4.64 ft 


5 ft/sec^ X X 6 ft/ft/sec^ =• 0.84 ft 

Total 4.5o RSS ed accelerometer errors is: 4.77 ft, or a 2a value of 

2.12 feet. 

C-1.2.2 GYRO ERRORS 

Errors existing prior to mode initiation: 

1 ft/ deg X 1.6® « 1.6 ft 

Error at align: 

(due to max. erection rats) 

1.0 ft 


a) Push over maneuver: 

b) Wind shears: 

Error Due to a^j term: 


RSS total 1.9 ft 


C-1.2.3 SUWIARY OF TOTAL ERRORS (Accelerometer and Gyro) 






for the Ky'« 2Ky case: 

Total accelerometer errors; 4.77 ft (4.5o) 
Total gyro errors; 1.9 ft (4. So) 

RSS total; 5.13 ft (4.5o) 

or, a total touchdown error of 2.28 ft on a 2o basis. 
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Figure C-2A. Rudder Limit and Crosswind Variations 
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Figure C-2B, Rudder Limit and Crosswind Variations 
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Figure C-4R No Rudder to Wheel Crossfeed 
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Figure C*SB. Doubled Ky and No Rudder to Wheel Crostfeed 
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Figure C-7B. Effect of MO OILS Discretization 
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Figure C-1 1. Effect of Wind Down Compentition tnd Inner RoU Loop G*in Vaiutions 
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Figure C-1 7. Response to 0.5 fps^ Error at h*1000 ft. 
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Figiire067. MODILS DiscretizatkMvif'TD C-81 
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Figure c*70. MODILS Discretization without Beam Noise, Y-rn 
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Figure 071. MODILS DiicretiMtion without Beam Noise , Y-nj 
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Figure 073. MODILS Discretization without Beaiin Noise. 4 jq 
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